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Preface 


This  study  was  prompted  by  a  phase  estimation  problem  of 
a  fading  signal  in  a  laser  line-scan  imagery  system.  The  sys¬ 
tem  is  sponsored  by  the  Air  Force  Avionics  Laboratory  and  it 
is  designed  to  produce  a  three  dimensional  image  of  the  ground 
surface  when  operated  from  an  aircraft. 

The  problem  of  nonlinear  phase  estimation  in  a  fading 
environment  is  quite  general  and  there  are  possible  applica¬ 
tions  in  many  areas  for  a  phase  estimator  which  is  insensitive 
to  amplitude  fluctuations.  The  wide  application  of  this  es¬ 
timation  problem  has  made  it  particularly  interesting  and 
stimulating  to  me. 

I  would  like  to  extend  my  thanks  to  my  thesis  advisor, 
Capt.  Stanley  Robinson , PhD,  MSEE,  BSEE,  for  his  guidance  and 
assistance  on  this  project  over  the  last  nine  months;  and  to 
my  committee  members  Professor  Peter  Maybeck,  PhD,  BSEE, 

and  Major  Gary  Reid,  PhD,  MS  ,  BS  for  their  very  useful 
comments  and  suggestions.  A  very  special  thank  you  is  also 
due  my  wife,  Maripat,  for  her  time,  encouragement,  suggestions, 
and  typing  throughout  this  project.  Her  help  has  meant  a 
great  deal  to  me. 
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Abstract 


An  airborne  laser  line-scan  imaging  system  can  determine 
the  relative  range  to  a  target  by  power  modulating  the  illum¬ 
inating  laser  and  detecting  the  relative  phase  of  the  reflect¬ 
ed  signal  with  respect  to  the  modulation  phase  reference. 

The  amplitude  of  the  reflected  signal  may  fluctuate  greatly 
due  to  such  factors  as  speckle  noise  and  the  varying  reflec¬ 
tivity  of  the  target,  resulting  in  the  severe  degradation  of 
performance  of  commonly  used  phase  detection  circuitry. 

In  this  report,  a  maximum-a  posteriori  (MAP)  phase  se¬ 
quence  estimator  is  developed  based  upon  the  Viterbi  algo¬ 
rithm.  The  signal  is  observed  in  additive,  zero  mean,  white 
Gaussian  noise  and  the  statistics  of  the  amplitude  are  as¬ 
sumed  to  be  Rayleigh.  Perfect  knowledge  of  each  realization 
of  the  signal  amplitude  is  assumed  and  the  sensitivity  to 
errors  in  the  amplitude  estimate  is  investigated.  Expres¬ 
sions  for  approximate  mean-squared  error  and  a  lower  bound 
on  mean-squared  error  are  developed  and  verified  by  a  Monte 
Carlo  simulation  of  the  estimator. 

For  quantization  of  the  phase  estimates  to  15  levels  and 
for  a  signal  with  constant  amplitude,  the  estimator  outper¬ 
forms  a  first  order  phase-locked  loop  (PLL)  by  2-6  dB  over 
the  range  of  input  signal  to  noise  ratios  of  3-17  dB. 

For  the  case  when  the  amplitude  is  rapidly  fading,  the 
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estimator  outperforms  a  first  order  PLL  by  greater  than  10  dB. 
This  performance  margin  is  based  on  the  Viterbi  estimator's 
performance  for  a  rapidly  fading  signal  compared  to  PLL 
performance  for  a  slowly  fading  signal.  Note  that  as  the 
amplitude  fading  bandwidth  of  the  PLL  signal  is  increased,  the 
PLL  will  lose  lock  completely,  so  that  the  .  r provement  margin 
presented  is  extremely  conservative. 

The  Viterbi  estimator  is  shown  to  be  insensitive  to  er¬ 
rors  in  the  estimate  of  the  signal  amplitude  for  input  sig¬ 
nal  to  noise  ratios  larger  than  2(3dB).  The  MAP  sequence  es¬ 
timator  is  also  shown  to  reduce  its  mean-squared  error  with 
increasing  decision  depth  (fixed  lag  of  the  Viterbi  algorithm), 
and  increasing  number  of  discrete  phase  levels.  The  performance 
depends  weakly  on  the  match  between  the  variance  of  the  actu¬ 
al  phase  and  the  variance  assumed  by  the  estimator. 


PHASE  SEQUENCE  ESTIMATION 
FOR  LASER  LINE-SCAN  IMAGERY 


IN  THE  PRESENCE  OF  RAYLEIGH  FADING 

I .  Introduction 

A  laser  line-scan  imagery  system  is  currently  being  de¬ 
veloped  by  the  Air  Force  Avionics  Laboratory  and  the  Environ¬ 
mental  Research  Institute  of  Michigan.  It  is  designed  for  air¬ 
borne  use  and  will  provide  a  three  dimensional  image  of  the 
earth's  surface.  Slant  range  information  is  obtained  by  mod¬ 
ulating  the  incident  laser  beam  with  a  periodic  signal  and 
then  comparing  the  phase  of  the  reflected  signal  with  a  known 
reference.  The  slant  range  to  the  surface  (assumed  to  be  flat 
earth)  can  be  combined  with  the  two  dimensional  reflected  sig¬ 
nal  power  image  to  provide  intelligence  data  to  a  human  oper¬ 
ator  or  to  automatic  cueing  equipment. 

The  signal  reflected  from  the  ground  target  is  perturbed 
in  amplitude  by  the  spatially  varying  reflectivity  of  the  tar¬ 
get  and  speckle  noise  (to  be  defined  subsequently)  (Ref6:4-5). 

This  amplitude  variation  detrimentally  affects  phase  estim¬ 
ation.  For  example,  a  linearized  first  order  phase-locked-loop 
will  lose  lock  as  the  fading  bandwidth  of  the  signal  is 
increased  (Ref  13). 

The  purpose  of  this  thesis  is  to  develop  a  phase  estimator 


which  is  insensitive  to  amplitude  variations.  The  performance 
of  the  estimator  will  be  calculated  and  results  will  be  verified 


by  a  simulation  of  the  estimator  algorithm. 

Background 

The  geometry  of  the  line-scan  imagery  system  is  depicted 
in  Figure  1.  The  aircraft  carrying  the  scanning  system  is 
shown  flying  in  the  Y  direction  while  the  laser  beam  is 
swept  across  the  ground  in  the  X  direction.  When  the  beam 
reaches  the  end  of  a  scan,  it  is  pointed  to  the  start  of  the 
next  scan  and  a  new  strip  of  ground  is  illuminated  as  the  beam 
scans  in  the  X  direction. 

The  beamwidth  of  the  laser,  speed  of  the  aircraft,  speed 
of  the  sweep,  and  altitude  of  the  aircraft  must  be  interrelated 
to  allow  the  scanner  to  sweep  contiguous  strips  of  ground  and 
cover  the  entire  area.  For  a  typical  beamwidth  of  one  milli- 
radian  (Ref3:2)  and  an  altitude  of  1000  feet,  the  illuminated 
area  is  a  spot  one  foot  in  diameter.  The  time  that  a  point  is 
illuminated  by  the  beam  is  the  dwell  time.  If  the  scan  length 
extends  to  1000  feet  on  either  side  of  the  ground  track  of  the 
aircraft,  and  if  the  aircraft  has  a  speed  of  600  knots,  the 
dwell  time  is  approximately  0.5  microseconds. 

A  block  diagram  of  the  scanner  system  is  shown  in  Figure 
2.  The  modulator  provides  a  reference  signal  to  the  phase  es¬ 
timator,  in  addition  to  modulating  the  power  of  the  laser 
source.  The  scan  optics  sweep  the  incident  beam  across  a 


Fig.  1.  Geometry  of  the  Laser  Line-Scan  System 

strip  of  ground  and  focus  the  reflected  signal  onto  the  op¬ 
tical  detector.  Either  direct  optical  detection  or  hetero¬ 
dyne  detection  is  suitable  for  the  receiver,  and  speckle 
noise  is  present  in  the  detected  signal  of  each  (Ref  6). 

Signal  processing  encompasses  any  amplification,  fre¬ 
quency  translation,  or  filtering  required  before  the  signal 
is  routed  to  the  phase  estimator  and  the  image  forming  cir¬ 
cuit  . 

Phase  Estimation  Problem 

As  described  previously,  amplitude  fluctuations  in  the 
received  signal  adversely  affect  the  task  of  phase  estimation. 
These  variations  in  amplitude  arise  from  three  main  sources. 


The  first  is  the  reflectivity  of  the  target  surface.  As 
the  illuminated  spot  scans  from  lightly  colored  objects  to  dark 
objects,  the  energy  in  the  reflected  signal  will  vary. 

The  second  cause  of  amplitude  variations  is  speckle  noise. 
When  light  is  reflected  from  a  surface  which  is  rough  with  re¬ 
spect  to  the  wavelength  of  the  incident  light,  a  random  ampli¬ 
tude  distortion  is  produced.  When  the  reflected  light  is  used 
for  image  formation,  the  distortion  manifests  itself  as  vary¬ 
ing  intensities  across  the  image.  This  is  called  speckle 
noise  (Ref  6:4-5). 

The  third  cause  of  phase  estimation  errors  is  the  receiv¬ 
er  hardware  itself.  The  signal  of  interest  has  a  large  dynam¬ 
ic  range  due  to  the  reflectivity  of  the  target  and  due  to 
speckle  noise.  Typically,  receivers  designed  to  accomodate  a 
large  dynamic  range  have  nonlinearities  in  their  phase  and  am¬ 
plitude  response  curves.  These  nonlinearities  can  distort  the 
phase  input  to  the  estimator  and  cause  additional  amplitude 
variations.  Also,  amplitude  modulation  (AM)  at  the  input  to  the 
receiver  may  be  converted  into  phase  modulation  (PM)  directly  by 
AM  to  PM  conversion. 

The  first  two  causes  of  estimation  problems,  reflectiv¬ 
ity  and  speckle,  are  inherent  in  the  nature  of  an  actively 
illuminated  system.  It  is  essential,  therefore,  to  use  a 
phase  estimation  scheme  which  is  robust  with  respect  to  am¬ 
plitude  fluctuations. 

Phase  estimation  problems  caused  by  imperfections  in  the 
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receiver  are  under  the  control  of  the  system  designer;  however 
an  estimator  which  is  insensitive  to  amplitude  variations  is 
still  desirable.  The  problem  of  a  receiver  design  with  good 
amplitude  and  phase  characteristics  can  be  treated  as  a  sep¬ 
arate  problem. 


Methods  of  Phase  Detection 


The  current  system  detects  the  zero  crossings  of  the  re¬ 
ceived  signal  and  the  reference.  The  relative  phase  is  deter¬ 
mined  by  measuring  the  time  between  zero  crossings.  This  type 
of  estimator  will  not  operate  below  a  threshold  set  by  the  par¬ 
ameters  of  the  circuit,  and  there  is  no  provision  for  making  a 
"reasonable"  estimate  during  periods  when  the  estimator  is  be¬ 
low  threshold. 

Phase-locked-loop  (PLL)  performance  is  well  known  for 
phase  estimation  applications  and  it  is  useful  as  a  benchmark 
against  which  other  estimators  can  be  measured.  A  PLL  exhibits  a 
threshold  effect  too,  usually  two  to  four  db  lower  than  that  of  a 
zero  crossing  phase  estimator.  Prior  statistics  of  the  phase 
process  may  be  available;  however,  PLL  estimators  do  not  usually 
make  use  of  this  prior  information. 

Both  zero  crossing  and  PLL  estimators  are  sensitive  to  ampli 
tude  fluctuations.  Amplitude  limiters  are  typically  used  to 
condition  the  input  signal  before  the  phase  is  estimated. 

The  phase  estimator  developed  in  this  thesis  is  one  which 
produces  the  maximum  a  posteriori  (MAP)  estimate  of  the  phase 
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sequence  after  some  fixed  delay.  The  algorithm  is  an  appli¬ 
cation  of  discrete  smoothing  and  employs  the  Viterbi  algo¬ 
rithm  (dynamic  programming)  for  actual  computation  of  the 
MAP  sequence. 

Scope  and  Assumptions 

In  this  thesis,  the  amplitude  of  the  received  signal  will 
be  considered  a  random  variable  with  Rayleigh  distribution  and 
independent  from  sample  to  sample.  The  amplitude  of  the  sig¬ 
nal  (before  addition  of  noise)  will  be  assumed  known.  Noise 
in  the  received  signal  will  be  assumed  to  be  white  and  gaussian. 
Phase  sequence  estimation  will  be  made  on  a  scan  by  scan  basis. 
The  phase  sequence  in  the  X  direction  of  Figure  1  will  be 
assumed  first  order  Gauss-Markov.  In  the  Y  direction,  the 
phases  are  assumed  independent . 

Overview 

This  report  begins  with  a  description  of  the  underlying 
system  constraints  and  the  nature  of  the  phase  estimation 
problem.  Chapter  II  develops  the  phase  model  and  algorithm 
for  a  MAP  phase  sequence  estimator  utilizing  the  Viterbi  al¬ 
gorithm.  The  performance  of  the  estimator  is  calculated  in 
Chapter  III  and  simulation  results  verifying  the  calculations 
are  given  in  Chapter  IV.  Comparisons  with  the  performance  of 
other  phase  estimators  are  made  in  Chapter  V  and  comments  on 
the  implementation  of  this  estimator  are  given  in  Chapter  VI. 
Finally,  conclusions  are  reviewed  and  recommendations  for  fur¬ 
ther  study  are  given  in  Chapter  VII. 
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II .  Theoretical  Background 

In  this  chapter,  the  theoretical  background  of  a  modulo-2ir 
MAP  phase  sequence  estimator  is  presented.  The  material  is  log¬ 
ically  divided  into  two  sections.  The  first  develops  a  model 
of  the  continuous  phase  process  to  be  estimated  and  the  discrete 
representation  of  that  process.  Statistics  of  the  process  are 
determined  and  an  expression  for  the  MAP  likelihood  of  a  sequence 
is  proposed.  Because  the  likelihood  measure  can  be  computed  re¬ 
cursively,  the  sequence  estimation  problem  can  be  solved  effi¬ 
ciently  through  the  use  of  the  Viterbi  algorithm  (Ref  9).  The 
application  of  the  algorithm  to  this  estimation  problem  is  dis¬ 
cussed  in  the  second  section  of  this  chapter. 

Signal  Model 

In  order  to  develop  a  model  of  the  phase  process,  it  is 
useful  to  begin  by  looking  at  the  signal  available  from  the 
direct  optical  power  detector. 

Optical  Detector  Output.  The  output  of  a  direct  optical 
detector  for  the  proposed  system  configuration  can  be  modeled 
by  (Ref  10:6) 

00 

y(t)  *  P0  j p(x)b(x-vt )m(t  -  ^ ) dx  (1) 

•00 

where  Pfl  is  the  peak  optical  power  including  fixed  optical 
losses,  p(x)  is  the  terrain  reflectivity  including  fluctu¬ 
ations  due  to  speckle,  b(x)  is  the  beam  filtering  function 
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(Ref  10:6),  m(t)  is  the  periodic  modulation  function,  r(x) 
is  the  slant  range  to  the  target,  v  is  the  sweep  velocity, 
and  c  is  the  speed  of  light.  Since  m(t)  is  periodic,  it 
can  be  expressed  as  a  Fourier  series  expansion: 


m(t)  =  yyyj 

n=-» 


Cne 


(2) 


where  f  is  the  fundamental  frequency  of  m(t)  and  the 
Cn's  are  the  Fourier  series  coefficients.  By  substituting 
equation  (2)  into  equation  (1)  and  taking  the  Fourier  trans¬ 
form,  the  frequency  domain  expression  for  Y(f)  can  be  ob¬ 
tained. 


Y(f>  -J'W t)>-  P.X)  B  [^-^]  Dn  [— O) 

n=-» 

In  equation  (3),  B'(f)  is  the  transform  of  b(x)  and  Dn(f) 

is  the  transform  of  p ( x ) expp f x— j  ,  the  terrain  informa¬ 
tion.  The  spectrum  of  the  terrain  information  is  shifted  to 
multiples  of  f  ,  the  fundamental  frequency  of  the  modulation 
signal,  as  shown  in  Figure  3. 

In  general,  the  lobes  of  the  spectra  can  overlap,  causing 
aliasing  errors;  however,  in  this  report  it  is  assumed  that 
there  is  no  aliasing  error.  It  should  be  noted  that  this  as¬ 
sumption  of  no  aliasing  error  usually  implies  that  filters 
with  "steep"  skirts  must  be  used  in  the  receiver.  For  many 
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Y(f) 


•Pis*  3.  Shifted  Spectra  of  Terrain  Information 

analog  filter  designs,  this  means  that  there  will  be  phase 
distortion  introduced  by  these  filters. 

If  the  signal  y(t)  is  passed  through  an  ideal  band¬ 
pass  filter  centered  at  nf^  ,  the  output  y'(t)  can  be 
described  as: 

00 

y'(t)  -  P0±g  Cn  J  p(x)b(x-vt)exp  £}27rnfm(t-  2r^x))] dx  (4) 
—00 

(Ref  10:7) 

The  filtered  signal,  y’(t)  ,  can  be  viewed  as  the  real 

part  of  a  compex  signal: 

y ’ (t )  =  Re[A(t)eJ<|,(t)e‘|2irnfnit]  (5) 
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The  complex  signal  can  be  related  to  equation  (4)  by 


y'(t)  =  A(t)cos[(j>(t)+27rnfmt] 

00  (6) 

=  P0  y* p(x)b(x-vt  )  I  Cj  2cos^n2irfm(t-  —  )+6nj  dx 

#  —00 


where 


C  = 
n 


C  I  e 
n  1 


je 


n 


(7) 


The  compex  signal  representation  will  be  useful  in  the  devel¬ 
opment  of  the  likelihood  ratio  derived  later  in  this  chapter. 
Without  loss  of  generality,  the  first  harmonic  of  fm  can  be 
selected  so  that  n=l 

Because  the  modulating  signal  m(t)  is  periodic,  the 
phase  <Kt)  can  be  measured  modulo-2u  with  respect  to  the 
reference  signal.  Since  the  slant  range  to  the  target  is  cal¬ 
culated  from  4> ( t )  ,  the  range  estimate  will  have  an  ambigu¬ 
ity  interval  of: 

ramb  "  27-  <Ref  3:19>  <8> 

m 


The  range  estimate  will  give  a  measure  of  the  range  within  the 
ambiguity  interval,  but  it  will  not  describe  in  which  ambigu¬ 
ity  interval  the  absolute  range  lies.  This  fundamental  limi¬ 
tation  allows  the  phase  to  be  estimated  modulo-2TT  with  no  loss 
of  range  information. 


Signal  Model.  The  signal  of  interest  will  be  observed  in 

additive  white  gaussian  noise  with  a  two  sided  power  spectral 

density  of  No/2  watts  per  Hertz.  The  task  is  to  estimate  a 

realization,  <j>(t),t  ft£tf  ,  of  the  phase  process  based  on  a 

s 

real  measurement 


j§(t)  =  ^(t)cos(jj>(t)+2irfmt)+J|I(t)  ts^t^tf+T  (9) 

where  jS(t)  is  the  real  measurement,  A(t)  is  the  amplitude 

of  the  signal,  JJJ(t)  is  white  gaussian  noise  and  the  "'v" 

denotes  a  random  variable.  The  interval  of  the  estimate  is 

from  t  to  t-  The  interval  of  the  measurement  is  allow- 
s  f 

ed  to  extend  past  tf  by  t  seconds.  By  allowing  the  measure¬ 
ment  interval  to  be  longer  than  the  interval  of  the  estimate, 
the  estimator  can  operate  with  a  fixed  lag  and  use  smoothing 
to  improve  the  estimate. 

For  the  scanning  system  described  in  chapter  I,  t 

s 

represents  the  start  time  of  a  scan  and  t^+x  is  the  time  at 

which  the  scan  ends.  In  practice,  an  estimate  of  £(t)  can 

be  made  for  the  entire  measurement  interval  t  <t£t-+x  :  how- 

s  i  ' 

ever,  the  estimate  of  ^(t)  for  t^St^t^+x  will  not  have  the 
benefit  of  the  full  lag  x 

The  continuous  time  model  of  equation  (9)  can  be  repre¬ 
sented  by  the  complex  discrete  time  measurement 


**  ^ke  +$k  k*l,2, . . . ,K+k( 


(10) 
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where  <j>k  =  <|>(t=kT)  on  the  interval  (k-l)TltlkT  ,  and  T 
is  the  sampling  interval. 

The  noise  term  in  equation  (11),  ,  is  a  white  Gaussian 

process  which  is  discrete  in  time.  It  is  assumed  that  £(t) 
is  constant  over  the  interval  (k-l)T5t5kT  and  £(t)=Ak  on 
this  interval.  The  amplitude,  £k  ,  is  assumed  to  be  a  ran¬ 
dom  variable  with  Rayleigh  statistics  and  independent  from 
sample  to  sample.  The  effects  of  a  correlation  between  ad¬ 
jacent  amplitude  samples  is  investigated  in  Chapter  IV.  The 
estimator  will  be  developed  with  the  assumption  that  the  se¬ 
quence  of  values  ,  k=l,2 , . . . ,K+Ko  ,  is  known.  The 

question  of  sensitivity  of  the  estimator  to  imperfect  know¬ 
ledge  of  £k  is  investigated  in  Chapter  IV. 

As  in  the  continuous  case,  an  estimate  of  the  sequence 

K  K+K  o 

{ 4>k )  i  will  be  made  based  on  the  measurement  sequence  {^k}j  , 

where  K0  is  the  fixed  lag.  In  the  discussion  of  the  Viterbi 

algorithm  later  in  this  chapter,  the  fixed  lag  Ko  will  be 

referred  to  as  the  decision  depth,  or  merely  depth. 

Phase  Statistics.  (Ref  9:8-11)  In  order  to  obtain  a 
stationary  distribution  on  the  process  ^(t)  ,  and  to  calcul¬ 

ate  conditional  transition  probabilities,  let  ^>(t)  be  a  ran¬ 
dom  walk  on  the  unit  circle,  where  ^(t)  can  take  on  values 
from  the  interval  (-tt,tt)  .  The  conditional  probability  of  a 

transition  to  <J>(t)  =  <j>.  given  that  <p(s)=<p  ,  t>s  ,  is 

t  s 

denoted  by  p(  <j> .  I  a)  .  This  conditional  density  satisfies 
the  partial  differential  equation  (Ref  9:8) 
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/ 


tr  P(»tl+S)  -  i0.a  Hr  PC+tUs)  (11) 

where  Oj  is  the  infinitesimal  variance  (random  walk  diffusion 

parameter)  and  -7 T<<p  ,<p  57r,t>s  .  The  boundary  conditions  are 

t  s 

lim  p(4>tUs)  =  6(^t— «*>s) 

P($t"-"|*s)  =  P(4>t=irUs)  (12) 


_9 

3*t 


p(<J>t=Tr  |  4>s) 


The  solution  for  equation  (11)  is  (Ref  9:9) 


p(»tUs> 


_ 1 

/2na0 2 (t-s) 


( 4>i--^.c;-n2Tr ) 2 1 
2 C7 0  2  (t-s)  ) 


(13) 


The  process  <|>(t)  is  thus  conditionally  approximately  gaussian 
with  mean  <P(s)  and  variance  afl2(t-s)  when  a20(t-s)  is 
small  (E{(4>t  -  <J>S ) 2  |  }  =  a2(t-s)). 

The  discrete  time  version  of  equation  (14)  can  be  obtain¬ 
ed  by  considering  the  sequence  C^}  arrived  at  by  sampling 
^(t)  at  intervals  t=kT  ,  k=0,l,2. . . .K.  The  transition  den¬ 
sity  from  to  is 


p(*ki*k.i>  - 


1 

/ZWoT* T 


^k"<*>k-l“27rn^2 

2a0zT 


(14) 
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/ 


Since  j£(t)  is  Markox,  {<}>k}  is  a  Markov  sequence  with  a 
probability  distribution  for  {<J>k}  x  of 

K 

PK({4>k}lK)  =  11  P^k^k-l5  (15) 

k=l 


The  initial  condition,  p(<j>il4>o)  ,  may  be  uniform  or  preset 

to  known  initial  condition. 

In  the  application  of  the  Viterbi  algorithm  to  this  extim- 
ation  problem,  the  phase  space  will  be  discretized  into 

values  £m  ,  m=0, 1, . . . .M-l  .  The  discrete  values  of  phase, 

£m  ,  will  be  selected  with  equal  spacing  from  the  interval 
(-tt,it)  and  M  will  be  restricted  to  an  odd  integer.  This 
allows  the  possibility  of  no  change  in  the  phase  from  the  cur¬ 
rent  time  to  the  next,  and  the  number  of  possible  levels  of 
phase  steps  in  the  positive  and  negative  directions  is  equal. 

It  is  necessary  to  calculate  the  transition  probabilities 
for  all  M2  pairs  j  and  m=0, 1, . . . ,M-1 

From  equation  (14)  let 


■  bp(*k'5mH>k-l=5j)  <16> 


where  b  is  chosen  so  that 

M-l 

X^k^irA-r^)  =  j=0,l,2,...M-l  (17) 

m=0 

The  sum  in  equation  (14)  must  be  truncated  at  some  level  which 
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gives  sufficient  accuracy.  If  the  discrete  phase  values  are 

equally  spaced  over  ( — tt  , tt  )  ,  then  the  value  p  depends  only 

on  the  distance  between  phase  values  |£.-£  I  modulo-27r  on 

1  j  m' 

( — rr , fr )  .  Therefore,  only  (M+l)/2  values  of  p  need  to  be 

stored  for  transition  probability  lookup,  as  opposed  to  the 
M2  number  of  possible  transition  pairs. 

Scharf  (Ref  9:12-13)  develops  the  relationship  between 
a  realization,  4>(t)  ,  of  the  continuous  phase  process  and 

the  discrete  time  realization  of  the  process  by  a  zero  order 
hold  (ZOH)  approximation,  4>k=<Kt=kT)  ,  (k-l)TltfkT  .  He 

shows  that  triples  (P{*},6,T)  can  be  selected  so  that  the 
ZOH  approximation  to  <f>(t)  is  close  to  4> ( t )  where  P{  • }  is 
the  probability  that  4>k  differs  in  value  from  <Kt)  by  more 
than  6  over  the  interval  T  .  Scharf  has  shown  that  this 
probability  is  given  by 

P{  *  >  =  4qI" — - — *1  (Ref  9:12)  (18) 

i/opTj 

where  Q[x]  is  the  Q-function 

00 

Q[x]  =  f  ~  exp(-52)d6  (19) 

J  /2rn 
x 

The  value  selected  for  T  depends  on  the  sampling  rep¬ 
resentation  of  Z^  and  the  speed  of  the  estimator  in  comput¬ 
ing  a  MAP  estimate.  In  any  event,  it  is  desired  that  T  be 
less  than  or  equal  to  the  dwell  time. 
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Phase  Sequence  Estimation.  The  real  and  imaginary  parts 
of  Zk=Uk+jVk  can  °^taine^  from  S(t)  by  a  quadrature  de¬ 
modulator  of  the  form  shown  in  Figure  4. 

This  is  equivalent  to 


Zk  = 


kT 

=  |  j  S(t)e 
(k-l)T 


- j2irf  t 
m 


dt 


(20) 


As  in  equation  (10),  Zk  is  the  sum  of  the  signal  ^Re 
and  white  gaussian  noise  with  a  power  spectral  density  of 
No/2  watts  per  Hertz.  If  and  <J>K  are  given  and  constant 

over  the  interval  (K-l)T^t^KT  ,  then  Uk  and  are  each 
conditionally  gaussian  with  zero  mean  and  variance  N0/T  (Ref 
9:15).  The  real  and  imaginary  parts  respectively  are  given  by 


«k-  AkcoS*k^ck 

(21) 

lkm  Aksl"»k+«sk 

where  jgck  and  ;jsk  are  the  discrete  time  narrowband  repre¬ 
sentations  of  the  noise  process  ^(t)  (Ref  14:502-504). 

A  likelihood  function,  ,  can  now  be  developed  in  order 

ir 

to  make  a  MAP  estimate  of  the  phase  sequence  {$k}j  The 

following  results  are  similar  to  those  obtained  by  Scharf  (Ref 
9:14-15)  but  have  been  expanded  to  include  the  case  where  the 
amplitude  of  the  signal,  Ak  ,  is  a  random  variable,  constant 
over  each  interval  T.  A  different  amplitude  value  is  allowed 
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Fig.  4.  Quadrature  Demodulator 


for  each  T  second  interval. 

The  likelihood  function  can  be  constructed 


D„  = 


K 

f o( 


(22) 


in  which  fi(*j*)  is  the  joint  conditional  probability  den- 

K  K 

sity  of  given  the  phase  sequence  and  fo(*|») 

|f 

is  the  joint  conditional  density  of  given  only  noise 


in  the  measurement. 


/ 


At  time  kT  ,  the  conditional  density  of  ^  given  (fi^ 


and  A. 


is 


k 


exp 


2  ' 

— 

zk'\e 

2N0 

"T" 

(23) 


and  the  conditional  density  of  z k  given  only  noise  is 


f(zk| noise)  exp 


.1! 


2Np 

T 


(24) 


Because  the  noise  is  white,  the  joint  conditional  densities 
fi  and  f0  can  be  obtained  by  appropriate  products  of  equa¬ 
tions  (23)  and  (24).  The  likelihood  ratio  (Refs  9:14  and  11: 
26-27)  portion  of  equation  (22)  (neglecting  the  priors  on  <j>) 
can  be  written 


/  T  \f  |  T  E 

fi_  \2TrN0/  exp|~  2N0  i=l 
JK  fo  /  m  \K  (  m  K 


Zk-Aie 


(srlr)7  expj-  sir  j 


(25) 


This  may  be  simplified  to 


or 


Lg  =  exp 


Lr  =  exp 


K 


2N0 


I 

i=l 


z.-A.e 


J*i 


K 

T  E 
i=l 


2N0  . 


K  (26) 

|“  5^7  S[_2ciAicos(,<;i"<l,i>+Ai2]|  <27> 
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Lk  =  exp 


/  K  K  x 

ZciAicos<n-*i>-  m  Z>4  '28> 

'  i=l  i=l  ' 


Because  the  IA^2  term  in  equation  (28)  is  independent  of  the 
data  of  interest,  it  will  be  neglected.  The  likelihood  func¬ 
tion  can  now  be  obtained  by  multiplying  by  the  priors  on 


{*k}  * 


(29) 


K  \  K 

dk  =  exp  ]h  2  C.A.cos^-^)  n  P(*i|*i_1) 

'  •  J  ^ _ -1 


The  expression  for  DK  gives  a  measure  of  the  likelihood 

V 

that  a  sequence  {<j>,  }  is  the  correct  sequence  given  the 

^  1 

K  K 

measurement  sequence  { Z.  }  and  the  amplitudes  {A.}  .  The 

K  i  1  1 

phase  sequence  which  maximizes  is  the  MAP  estimate  of 

‘*k>?  • 

By  taking  the  natural  logarithm  of  DK  ,  this  result  can 

be  simplified  further  to 

K  K 

rK  =  ln(DK)  =  j^]>]AiCiCOs(V<})i)  +  [p(*i  I  (30> 


where  Ci  is  the  magnitude  of  the  measurement  Z±  and 
is  the  argument  of  Zi  ,  as  shown  in  Figure  4.  The  subscript 
K  on  the  likelihood  measure  rK  indicates  that  rK  has  been 

V 

calculated  for  the  sequence  of  measurements  {Z^j  .  In 
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general,  a  I\  may  be  calculated  for  measurements  from  Zj 
to 

A  more  useful  expression  can  be  obtained  by  noting  that 
Ti  can  be  computed  recursively  and  that  the  phase  process  is 
first  order  Markov. 


r. 

i 


T.  A  .C.  cos(ij) .  — 4> .  )+ln 

l-l  No  11  vri  ri' 


P(4>i  1 4>i__i ) 


(31) 


The  recursive  expression  for  leads  very  nicely  to  the 

application  of  the  Viterbi  algorithm  for  phase  sequence  estim¬ 
ation. 

Viterbi  Algorithm 

As  previously  stated,  the  desired  result  is  selection  of 

A 

the  sequence  out  of  all  possible  sequences  which  max- 

K 

imizes  given  the  measurement  sequence  (Z^)^  .  This  is 

the  MAP  estimate  of  the  phase  sequence.  If,  at  each  time  k  , 

the  phase  4>k  can  take  on  one  of  M  values,  then  there  are 
K 

M  possible  sequences  to  search  through  in  order  to  find  the 
one  that  maximizes  .  The  complexity  of  this  calculation 

quickly  becomes  too  large  to  be  useful  as  M  and  K  become 
large. 

Because  rR  can  be  computed  recursively,  as  in  equation 
(31),  this  exponential  increase  in  complexity  can  be  avoided. 
It  will  be  convenient  to  define  as  the  path  metric 

VJ',n>  “  It  Akckcos<tk-'V5J)+li>  PCVCjM’k-l-V  (32) 
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therefore 


rK  =  rK-l+YK(:i  ,m) 


(33) 


A  £ 

The  sequence  1<J>  ^  }  x  which  maximizes  r  K  while  passing 

A  A 

through  $  ^  on  its  way  to  cf>K  must  do  so  by  following 

A  K— 2 

sequence  {$k}#~  which  maximizes  rk-1  .  If  it  did  not,  a 

a  K-2 

better  sequence  ($k}i""  could  be  chosen.  Therefore,  at  time 

f  A  V 

K  ,  it  is  only  necessary  to  store  M  sequences  , 

one  ending  in  each  of  the  M  phase  values  at  time  K  ;  M 
values  of  rK  corresponding  to  the  M  sequences;  the  current 
time  counter,  K  ;  and  the  transition  probability  values  (Ref 
9:24-27). 

To  extend  the  MAP  estimate  of  phase  to  time  K+l  ,  the 
estimator  algorithm  must  take  the  new  measurement  ZK+1  and, 
for  each  phase  value  *  j=0,l,...M-l  ,  calculate 

rK+i  for  each  sequence  plus  <J)K+i=5j  •  For  each 

value  of  j  ,  the  maximum  rK+^  is  selected  and  the  value  of 

t r 

JK+1=5j  is  concatenated  with  the  sequence { <f>k )i  which  ends  in 
$k=5m  '  At  the  enci  of  the  calculation,  zk+1  wiH  no  longer 

be  required  and  it  can  be  discarded.  There  will  be  M  phase 
sequences  stored,  one  ending  in  each  of  the  M  possible  dis¬ 
crete  values  of  phase.  The  sequences  can  be  extended  another 
step  in  time  by  taking  measurement,  zj^+2  »  and  repeating  the 
calculations. 

This  maximization  algorithm  can  be  stated 
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/ 


(Ref5 :24)  (34) 

or  in  algorithm  form 

For  each  time  k=l,....K 
Begin 

Obtain  the  current  measurement  Zk 
For  each  phase  value  ?.  ,  j=0,l,...M-l 

J 

Begin 

For  each  phase  value  ,  m=0,l,...M-l 

Begin 

Calculate  rk  ssrk_1+Yjc(  j  ,m) 
m 

End 

Pick  the  biggest  rk 
/v  m 

Concatenate  with  the  sequence  ending 

in  *k-r5™ 

^  K 

Set  rk  =  rk  for  the  sequence  { <t> k } i  ending 

m 

in  Cj 
End 
End 

At  each  time  k  ,  there  will  be  a  sequence  of  phase  val¬ 
ues  ending  in  each  of  the  M  possible  phase  values.  These 
are  called  survivor  sequences.  In  general,  there  will  be  sur¬ 
vivor  sequences  from  the  previous  time  (k-1)  which  are  not 


Max 

{$} 


K 

K-1 


“iXK-2 [rK-1+YK( J ’m)] 


<*), 
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concatenated  with  a  phase  value  to  form  a  current  survivor 
sequence.  These  unused  sequences  need  not  be  retained. 


The  phase  values  stored  in  the  M  sequences  tend  to  con¬ 
verge  to  a  single  value  if  the  sequences  are  inspected  far 
enough  back  in  time.  Because  of  this  convergence,  at  time 

K  it  is  possible  to  determine  the  MAP  estimate  of  the  phase, 

/V  A  K 

4>K_k0  >  by  selecting  the  sequence  with  the  largest 

A 

IV  and  outputting  the  <J)V  value  from  that  sequence.  If 

xi. — iv  o 

all  the  M  sequences  have  not  converged  to  a  single  value 
of  <|>K  ,  then  the  ambiguity  can  be  resolved  by  making 

K0  larger.  It  should  be  noted  that  this  does  not  guarantee 


that  is  tbe  closest  neighbor  to  the  actual  phase 

A 

<j>K_K  ^  ,  only  that  is  the  MAP  estimate  of  the  true 

value . 


The  number  K0  is  the  fixed  lag  of  the  system  and  is  re¬ 
ferred  to  as  depth,  or  decision  depth. 

The  complexity  of  the  estimator  can  be  seen  from  the  al¬ 
gorithm  listing.  If  the  times  required  for  concatenation, 
picking  the  biggest  rK  ,  and  updating  the  new  values  of  r 
are  ignored,  then  the  majority  of  time  is  spent  in  calculation 

of  the  values  of  (r^  is  the  mth  value  computed  in 

m  m 

the  Viterbi  Algorithm).  They  must  be  calculated  M  times  for 
each  of  the  M  values  of  phase,  at  each  time  k=l,2,...K. 
Therefore  the  complexity  of  computation  is  proportional  to 

tr 

KM2  .  This  is  better  than  the  M  complexity  required  by 
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the  brute  force  global  search  method  mentioned  at  the  beginning 
of  this  section.  In  addition,  the  Viterbi  algorithm  estimator 
outputs  phase  value  estimates  with  a  fixed  lag  in  time,  as  op¬ 
posed  to  the  global  search  method  which  must  store  the  entire 
measurement  sequence  and  then  calculate  the  entire  MAP  estim¬ 
ate  of  the  phase  sequence. 

A  signal  model  has  been  developed  which  describes  the  re¬ 
ceived  signal  of  the  laser  line  scan  system.  The  phase  of  the 
received  signal  can  be  represented  by  a  sequence  of  values 
over  time.  Using  the  Viterbi  algorithm,  it  is  possible  to 
make  a  maximum  a  posteriori  estimate  of  this  phase  sequence. 

The  question  of  how  well  this  estimator  performs  is  addressed 
in  Chapter  III. 
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Ill  Performance  Analysis 


An  analysis  and  prediction  of  the  performance  of  an  est¬ 
imation  scheme  is  necessary  in  order  to  determine  if  the  est¬ 
imator  is  worth  building  and  to  compare  its  performance  with 
competitive  estimators.  In  this  chapter,  the  mean-squared 
error  (MSE)  of  the  estimate  is  presented  as  the  measure  of 
quality  for  the  MAP  phase  sequence  estimator.  By  using  the 
concept  of  error  events  and  minimum  free  distance,  approximate 
values  of  MSE  are  predicted  for  various  signal  energy  to  noise 
power  density  ratios  and  number  of  levels  of  phase  discretiza¬ 
tion.  Also,  a  strict  lower  bound  on  the  MSE  is  calculated. 

The  question  of  bias  in  the  estimate  is  discussed  and  the 
limitations  of  the  predicted  MSE  are  presented. 

Mean-Squared  Error 

A  V 

The  sequence  of  phase  values  {$^1  which  approximates 
the  true  phase  <p(t)  ,  0<t£KT,  can  be  viewed  as  a  path 

through  a  trellis.  Figure  5  shows  one  such  path  through  a 
trellis  in  which  there  are  five  equally  spaced  discrete  phase 
values.  The  value  of  <p(t)  ,  (K-l)T<t£KT  ,  is  represented 

by  the  node  at  K  through  which  the  path  passes.  The  lines 
connecting  nodes  represent  transitions  from  one  phase  value 
to  another. 

In  this  example,  the  phase  is  shown  changing  linearly 
from  each  point  in  time  to  the  next.  The  linear  phase 
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Fig.  5.  Phase  Trellis,  M=5 


transition  case  and  a  step  phase  transition  case  will  both  be 
used  in  calculating  the  predicted  mean-square  error  of  the 
estimator. 

The  first  case  can  be  viewed  as  a  piecewise  lineariza¬ 
tion  of  the  continuous  phase  (range).  This  is  a  reasonable 
(although  simplified)  model  for  the  modulo-2ir  phase  which 
is  proportional  to  the  slant  range  to  the  terrain  being 
scanned.  This  linear  transition  case  will  be  useful  in  pro¬ 
viding  a  prediction  of  the  mean-squared  error  (MSE)  of  the 
system  when  the  input  phase  is  a  result  of  a  physical  pro¬ 
cess,  such  as  scanning  of  terrain. 
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The  second  case,  step  phase  transition,  is  included  be¬ 
cause  this  case  matches  the  model  for  the  Monte  Carlo  simula¬ 
tion  performed  in  Chapter  IV.  In  the  simulation,  the  input 
phase  sequence  is  generated  as  discrete  values  which  are  con¬ 
stant  over  a  T  second  interval.  This  corresponds  to  a 
step  change  in  phase  in  each  interval. 

The  free  distance  (to  be  defined  later)  between  two  sig¬ 
nals  in  signal  space  is  smaller  for  the  linear  transition 
case  than  for  the  step  transition  case.  It  will  be  seen,  in 
the  following  sections,  that  this  results  in  a  higher  prob¬ 
ability  of  making  an  error  in  phase  estimation  for  the  linear 
transition  case. 

In  this  report,  the  correct  path  is  defined  as  the  se¬ 
quence  of  phase  values,  which  is  closest  in  absolute  value  at 
each  discrete  time  kT  to  the  actual  phase  value.  Conversely, 
an  error  path,  ei  ,  is  the  i’th  path  which  is  not  identical 
to  the  correct  path  at  all  k«(l,2, . . ,K).  For  a  sequence  of 
length  K  which  can  take  on  one  of  M  discrete  values  at  each 
time  k  ,  there  will  be  one  correct  path  and  M  -1  error 
paths.  The  entire  collection  of  error  paths  is  denoted  as  E  . 
It  will  also  be  useful  to  define  P£  as  the  probability  that 
the  estimator  has  chosen  the  correct  path,  that  is,  the  se¬ 
quence  of  values  output  from  the  estimator  is  the  same  as  the 
sequence  of  values  of  the  correct  path.  In  a  similar  manner, 

P  is  defined  as  the  probability  that  the  estimator  has  chosen 
ei 
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error  path  e.  and  P_  is  defined  as  the  probability  that  the 

*  l  hi 

estimator  is  on  any  error  path. 

It  follows  that 


(35) 


and 


(36) 


The  estimator  chooses  either  the  correct  path  or  an  error 

path.  These  two  events  are  mutually  exclusive  and  exhaustive, 

< 

therefore  the  mean-squared  error  (MSE)  at  any  point  in  time  on 

l 

the  trellis  can  be  defined  as 

MSE  =  (MSE|C)Pc+(MSE|E)Pe 

-  (MSE|C)(1-Pe)+(MSE|E)Pe  (37) 

=  (MSE|C)(1-Pe)+  £  (MSE|e.)Pe 

eieE  1 

where  (MSE|C)  is  the  mean-squared  error  given  the  estimator 
is  on  the  correct  path. 

The  ( MSE | C )  term  is  due  to  the  quantization  noise  of  the 
approximation  of  4>(t)  by  one  of  a  set  of  M  discrete  levels 
on  the  interval  (-7r,Tr)  .  If  it  assumed  that  <f>(t)  is  dis¬ 
tributed  uniformly  over  each  quantization  interval,  then 

( MSE  |  C )  =  —  (38) 

3M2 
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The  assumption  of  uniformity  over  a  quantization  interval  is 
reasonable  if  M  is  sufficiently  large  and  the  probability 
density  function  (pdf)  of  the  phase  is  "smooth"  so  that  the 
pdf  of  ( t )  is  essentially  constant  over  each  quantization 
interval . 

The  (MSE|E)  is  similarly  defined  as  the  mean-squared 
error  given  that  the  estimator  has  chosen  any  path  other  than 
the  correct  path.  The  average  mean-squared  error  at  a  point 
in  time  (given  that  the  estimator  has  chosen  error  path  e^) 
is  denoted  (MSE|e^)  .  The(MSE|e^)  is  averaged  over  the 
length  of  the  error  event.  It  will  be  seen  in  the  following 
development  that,  for  many  cases  of  interest,  (mse|e)  is  dom¬ 
inated  by  one  (MSEje^J  term  in  equation  (37).  The  dominant 
term  has  a  length  of  one  and  the  MSE  defined  in  equation  (37) 
becomes  a  point  MSE  for  these  cases. 

For  the  calculation  of  (MSE|e^)  it  will  be  useful  to 
use  the  following  definitions  for  the  signal  of  interest  and 
its  phase  (Ref  1:42.4.1) 

s(t,a)  =  A  cos(27rfmt  +  ({>(a) )  (39) 

co  n 

4>(s0  =  J  aig(a-iT)da-«<t<(n+l)T  (40) 

—00  i  =—00 

where  M  is  the  number  of  discrete  phase  quantization  levels, 
and  a  is  a  sequence  of  integers  selected  from  the  alpha¬ 

bet  (0,±l,±2,...±(^i)  ). 
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For  convenience  in  the  development  of  the  free  distance  later 

\ 

in  this  section,  the  amplitude  is  defined 

>=VF  («> 

where  E  is  the  energy  in  the  signal  s(t,a)  over  a  T  sec¬ 
ond  interval. 

The  function  g(t)  is  selected  so  that 

00 

/  g(t)dt  =  1  (42) 

—  00 

The  shape  of  g(t)  determines  how  the  phase  of  s(t,a)  pro¬ 
gresses  from  one  node  on  the  trellis  to  the  next.  If  g(t)=i  , 
0<t<T  and  g(t)=0  elsewhere,  then  the  phase  changes  linearly 

with  time.  If  g(t)=<5(t)  ,  then  the  phase  makes  a  step  jump 

2ira . 

of  height  — jjj —  in  each  T  second  interval.  The  elements  ou 

2  7T 

of  a  specify  the  number  of  jp  radian  steps  in  phase  the  sig¬ 
nal  changes  on  the  interval  iT<t5(i+l)T,  The  ou  are  limited 

M-l 

in  magnitude  to  — to  reflect  the  modulo-2iT  nature  of  the 
estimator. 

Selection  of  a  function  for  g(t)  which  gives  a  linear 
phase  change  (or  similar  smooth  change  in  phase  versus  time) 
is  useful  for  modeling  a  physical  svstem  such  as  the  laser 
line  scanner  described  in  the  Introduction.  This  is  due  to 
the  bandwidth  constraints  which  limit  the  dynamics  of  the 
phase  change.  The  simulation  results  presented  in  Chapter  IV 
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are  based  on  an  estimator  model  in  which  the  measurements  are 
discrete  in  time  and  quantized  in  value.  Thus,  the  Dirac 
delta  function  more  closely  describes  the  phase  transitions 
for  the  estimator  simulation. 

For  the  rest  of  this  section,  a  is  a  sequence  of  in¬ 
tegers  which  describes  the  correct  path  of  the  estimator  and 
the  sequence  £  will  represent  one  of  the  error  paths  which 
differs  from  a  at  L  nodes.  The  value  L  is  defined  as 
the  length  of  the  error  event.  An  example  of  a  correct  path 
a  and  an  error  event  £  with  L=2  is  shown  in  Figure  6. 

In  this  example  of  an  error  event  the  error  path  diverges 
from  and  merges  with  the  correct  path  only  once.  In  general, 

£  could  diverge  from  the  correct  path  several  times  and  it 
need  not  merge;  however,  it  will  be  seen  in  the  following 
development  that  it  is  the  L=1  "minimum  error  path"  which 
is  of  primary  interest. 

The  averaged  MSE  given  error  path  e^  (described  by  £  ) 

and  correct  path  a  can  be  defined  as 

A  -i 

(MSE|ei)  =  £  2J  [<Kaj)-<KB;j)]2  (43) 

J 

where  the  sum  is  taken  over  the  L  nodes  at  which  a  differs 
from  £ 

Forney  has  shown  (Ref  2:273)  that  error  events  are  recurrent, 
or  statistically  independent  of  each  other.  The  averaged  MSE 
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Fig.  6.  Error  Event 

for  a  particular  error  event  is  defined  as  the  average  of  the 
point  mean-square  errors  over  the  length,  L  ,  of  the  error 
event.  It  will  be  seen  that  the  most  probable  error  path 
has  length  L=1  ,  and  that  (MSE|ei)  for  this  most  probable 
error  path  collapses  to  a  point  mean-square  error. 

If  the  error  path  diverges  once  at  time  (n-L-l)T  and 


merges  once  at  time  nT  , 

then 

(MSE)e^)  =  (|I)!  i 

n-1 

E 

j=n-L 

*  J 

E  ai- 

i=n-L 

3  i* 

£  »i 
i=n-L 

(44) 
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Further  development  of  the  total  mean-squared  error  re¬ 
quires  an  expression  for  the  probability  of  an  error  path. 

Probability  of  Error  Path 

The  desired  error  measure  is  the  probability  per  unit 
time  that  the  estimator  will  select  a  particular  incorrect 
node  in  the  phase  trellis.  From  this,  the  MSE  per  unit  time 
could  be  calculated  using  an  equation  of  the  same  form  as 
equation  (37).  This  probability  measure  is  not  directly 
available;  however,  useful  approximations  and  bounds  to  the 
mean-squared  error  can  be  obtained  by  calculating  the  prob¬ 
ability  that  an  error  event  starts  at  time  nT  given  that 
the  estimator  was  on  the  correct  path  at  time  (n-l)T  (Ref  4 
273-274). 

The  probability  that  the  estimator  selects  a  particular 
error  path  instead  of  the  correct  path  is  a  binary  hypothesis 
problem  and  depends  on  the  euclidean  distance  between  the  two 
signals  in  signal  space.  These  error  events  are  not  mutually 
exclusive,  therefore,  the  probability  that  any  error  event 
occurs,  P£  ,  can  be  upper  bounded  by  the  sum  of  all  the  bi¬ 
nary  decision  error  probabilities  (union  bound).  A  lower 
bound  on  the  probability  of  any  error  event  is  the  largest  of 
the  probabilities  of  specific  error  events  (Ref  4:273-274). 
These  bounds  are  summarized  as  follows  (Ref  4:274) 

T  N£pe£?pn -  MrN+£%  (45) 
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The  binary  hypothesis  probability  of  an  error, 
(Refl:42.4.2) 


P 


e 


i 


=  Q 


££ jki&y 
/5n7 


,  is 


(46) 


where  N0  is  the  one  sided  power  spectral  density  for  additive 
white  gaussian  noise,  D  is  the  euclidean  distance  between  two 
signals  defined  by  a  (the  correct  path)  and  j3  (the  i'th 
error  path),  and  Q( • )  is  the  error  function  defined  by 

2 

2 ?  -y 

Q(x)  =  -i-  /  e  2  dy  (47) 

/2tt  X 


The  Q  function  is  monotonically  decreasing  and  rapidly 

becomes  small  as  D  becomes  large.  As  a  result,  the  P  in 

E 

equation  (45)  can  often  be  approximated  by 


P 


E 


where  D  is 
m 


D 


m 


(48) 


(49) 


Note  that  Pj/-PE  by  equation  (45).  The  minimum  is  taken 
over  all  possible  error  paths  e^ 

The  total  mean- square  error  can  now  be  approximated  by 

MSE  "  W  (1-pP+<MSElem>Pi  (50) 
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In  equation  (50),  em  is  the  error  path  which  is  the  minimum 
distance  from  the  correct  path  defined  by  a 

Since  is  a  lower  bound  on  P£  ,  a  lower  bound  on  MSE 

can  be  defined  as 

MSE>(MSE|em)P"  (51) 

An  upper  bound  on  the  MSE  is  desirable;  however,  this 
would  require  evaluation  of  the  union  bound  sum  in  equation  (45). 
This  would  provide  an  upper  bound  on  P£  and,  therefore,  an 
upper  bound  on  MSE  by  a  substitution  into  equation  (50). 

The  large  number  of  paths  over  which  the  sum  must  be  cal¬ 
culated  makes  evaluation  of  an  upper  bound  infeasible,  so  an 
upper  bound  on  MSE  is  not  evaluated  in  this  report.  As  will 
be  shown  in  the  following  development,  a  measure  of  the  "tight¬ 
ness”  of  the  approximation  of  MSE  is  the  value  of  the  argument 
of  the  Q  function  used  in  equation  (48).  This  will  give  an 
indication  of  the  quality  of  the  approximation  without  actual 
calculation  of  an  upper  bound. 

The  minimum  free  distance,  ,  must  be  calculated  be¬ 

fore  actual  values  of  MSE  can  be  evaluated. 

Minimum  Free  Distance 

In  this  section,  the  normalized  free  distance,  d=D/y/21? 

(Ref  11:42.3.2),  will  be  developed  for  two  cases  of  g(t) 
and  then  the  minimum  free  distance,  D  ,  wil2  be  evaluated. 

The  signal  and  phase  are  modeled  as  in  equations  (39)  and  (40), 


y 


respectively.  It  should  be  noted  that  the  amplitude  of  the 
signal  is  constant  in  this  development.  The  case  of  varying 
amplitude  will  be  discussed  in  a  following  section. 

Linear  Phase  Transition.  For  the  first  case  of  interest, 
the  phase  changes  linearly  from  one  time  interval  to  the  next. 
This  case  can  be  viewed  as  a  piecewise  linearization  of  a  con¬ 
tinuously  changing  phase. 

The  phase  4>(°0  is  zero  for  t<0  and 

g(t)  =  i  0<t<T  (52) 

-  0  elsewhere 
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the  phases  differ.  If  oi  T>>1  .  then 

m 


(j+l)T 


£  / 

j-n-N+1 


M  S  (ai“3i)+  T(aj+l“Sj+l)  " 
i=l 


(55) 


dt 


The  integral  in  equation  (55)  must  be  evaluated  in  two 
cases.  When  a  ...  =  8  .  , 

j  +  i  j+i 


y**dt  =  T  COS  yip  £ 


(56) 


When  aJ+1*0J+1 


dt  = 


2ir(aj+l-Bj+l) 


sin  I1  E  (ai-ei)+(aj+i-ej+i) 


(57) 


j  I 

sin  r  E(0'i'ei) 

i=l 

The  term  J- dt  represents  the  entire  integral  in  equation  (55). 

From  equations(55) ,  (56)  and  (57),  it  can  be  seen  that  an 
error  path  defined  by  jj  will  have  a  minimum  distance  from  the 
correct  sequence  when  its  phase  differs  from  a  for  only  two  T 
second  intervals.  This  is  illustrated  in  Figure  7,  with  N=2 


/ 


4tt 

5“  * 


Fig.  7.  Minimal  Error  Path  _0 

and  the  phase  difference  occurring  at  the  beginning  of  the 
sequence.  The  location  of  the  phase  difference  in  time  does 
not  affect  the  value  of  d2 
For  this  error  event, 

ai+a2  =  ®l+^2  (58) 

therefore 

a  -B  =  -(a  -B  )  (59) 

11  2  2 

so  the  resulting  normalized  squared  distance  is 

d2  =  2  jl-sinc  jp  (ai-Bi)l  (60) 
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(61) 


where 

sino(x)  - 

X 

To  minimize  d2  in  equation  (60),  the  argument  of  the 
sine  function  must  be  made  as  small  as  possible.  For  a  fixed 
value  of  M  ,  this  can  be  done  by  selecting  .  (They 

cannot  be  equal  because  that  would  imply  that  a=j3  ,  and  no 
phase  difference  existed.) 

The  minimum  free  distance  squared  can  now  be  stated  as 

D2  =  2Ed2  =  4E(  1-sinc  §2-)  (62) 

m  m  '  M  '  v  ' 


It  should  be  noted  that  for  each  correct  path  defined  by 
a  ,  there  are  two  error  paths  which  have  the  same  minimum  dis¬ 
tance,  Dm  ,  from  the  correct  path.  Equation  (48)  is  modified 
to  include  both  of  these  minimum  distance  paths. 


P 


(63) 


This  result  can  be  combined  with  equations  (44) , (50) , (51) 

and  (62)  to  calculate  the  approximate  MSE  and  a  lower  bound  on 

the  MSE  as  shown  in  Figures  (8)  and  (9).  The  MSE  values  in  all 

2 

figures  have  been  divided  by  to  normalize  them  with  re¬ 
spect  to  a  uniform  phase  distribution  over  ( — tt , tt )  :  MSEN  * 

(MSE)  |r  . 
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Figure  8  Predicted  MSEN,  Linear  Transition,  Constant 


Figure  9  Lower  Bound  on  MSEN,  Linear  Transition,  Constant 


The  equation  for  approximate  MSEN  is 


MSEN 


(64) 


The  lower  bound  is  initially  calculated  by 


(65) 


Since  the  MSEN  cannot  be  less  than  the  quantization  noise  due 
to  M  levels,  the  lower  bound  depicted  in  Figure(9)  is  actual¬ 
ly  calculated  by 


MSE 


24  ~ 
»  °r>  M2  ^ 


(V2i7  (1-sino  sr) ) 


(66) 


Figures  (10)  and  (11)  display  the  predicted  MSEN  and  low¬ 
er  bound  on  MSEN  versus  an  abscissa  of  E/No  .  As  shown, 
the  performance  of  the  estimator  is  limited  by  the  quantization 
noise  caused  by  representing  a  continuous  phase  by  one  of  M 
discrete  values.  The  flat  portions  of  these  curves  denote  the 
domain  of  E/No  over  which  quantization  noise  is  the  limiting 
factor. 

Step  Phase  Transition.  The  second  case  of  interest  is  that 
in  which  the  phase  is  allowed  to  make  a  step  change  at  time 
t=kT  ,  for  k=0, 1,2, ... ,K  .  In  this  case,  g(t)=  6(t)  ,  the 

Dirac  delta  function.  As  previously  described,  this  case  models 
the  Monte  Carlo  simulation  described  in  Chapter  IV. 
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Figure  10  Predicted  RISEN,  Linear  Transition,  Constant 


LINEAR  PHASE 


1 


Lower  Bound  on  MSEN,  Linear  Transition,  Constant 


The  phase  function  is 


J. 


The  lower  bound  cannot  be  less  than  the  quantization  noise, 
therefore 


MSE  >  Max 


1 

M2 


24  n 

or.  m  « 


(72) 


Figures  (12)  and  (13)  show  the  results  of  the  Normalized 
MSE  calculations  versus  M(quantization  levels)  for  the  step 
transition  case.  The  predicted  MSEN  results  for  this  case  are 
slightly  better  than  those  predicted  for  the  linear  phase 
transition  case.  Equations  (70)  and  (72)  are  also  used  to 
calculate  MSEN  and  a  lower  bound  on  MSEN  versus  E/Nq  as 
displayed  in  Figures  (14)  and  (15). 

As  in  the  linear  transition  case,  the  flat  portions  of 
the  curves  in  Figures  (14)  and  (15)  indicate  regions  where 
the  quantization  noise  dominates  the  MSE  performance  of  the 
estimator. 


Random  Variable  Amplitude 

The  predicted  mean-square  error  shown  in  Figures  (8) 

through  (15)  are  all  based  on  a  constant  amplitude  signal.  If 

2E 

the  amplitude,  £=  _Z  ,  is  a  random  variable  then  the  signal 
can  be  defined  as 


s(t 


.=>  -V? 


cos(2Trfmt+<j)(t ,  a)) 


(73) 


It  is  assumed  that  ^  is  constant  over  each  T  second 
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Figure  12  Predicted  MSEN,  Step  Transition,  Constant 
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Figure  13  Lower  Bound  on  MSEN,  Step  Transition,  Constant 


Figure  14  Predicted  MSEN,  Step  Transition,  Constant 


STEP 


Figure  15  Lower  Bound  on  MSEN,  Step  Transition,  Constant 


/ 


interval  and  E^=JS(t)  for  (i-l)T<t<iT  .  Perfect  knowledge 
of  the  Ei's  by  the  estimator  is  assumed  and  the  E^ 1 s  are 
assumed  identically  independently  distributed.  The  signal 
can  now  be  represented  by. 

jsr 

s(t,ai,Ei)  =  cos(2irfmt+(j>(t,ai))  (74) 


(i-l)TltliT 


A  signal  with  phase  sequence  j5  is  defined  to  have  .amplitude 
values  B.^ 

Calculation  of  a  predicted  MSE  for  the  varying  amplitude 
case  requires  evaluation  of  the  minimum  euclidean  distance  Dm 
for  the  phase  transition  being  considered  and 


(75) 


The  notation  E{*}  is  the  expected  value  operator  on  the  ran¬ 
dom  variable  E  and  D  is  a  function  of  the  random  variable 

^  ^m  j/ 

J|g  .  Due  to  the  complexity  of  calculating  Ej 
closed  form,  the  Q [ • ]  function  is  approximated  by 


e  random  van; 

■■ 


Q[x]  t  be~2 


(Ref  11:40)  (76) 


This  allows  calculation  of  the  expected  value  of  Pg" 
which  is  now  defined  as 
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E{P-)  -  e  j*  expfJr- 


(77) 


The  predicted  mean-squared  error  can  now  be  calculated 
for  the  two  phase  transition  cases  previously  described. 

Linear  Phase  Transition.  The  normalized  squared  dis¬ 
tance  is  defined  as  in  equation  (54);  and  g(t)  is  for 
O^t^T  ,  and  zero  otherwise. 

For  the  error  event  when  N=2  ,  a  1+012  =  61  +  32  and 

;  it  can  be  shown  that 


D2  =  2(]g1+jg2)-2jgxsinc(g^al-B1))- 


(78) 


2E2sinc(|^(a1-81)) 


The  minimum  free  distance  occurs  when  |ai-8il=l  ,  as  in 
the  previous  examples. 

The  amplitude  of  the  signal,  £  ,  is  assumed  to  have  a 

Rayleigh  probability  distribution  function.  A  transformation 
of  variables  (Ref  7:129-131)  gives  the  probability  density 
function  for 


-2aE 

fjg(E)  =  I2-  e  T  ,  E>0 


(79) 


,  otherwise 


where 


E(A2}  =  - 


(80) 
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The  expected  value  of  P£  over  and  is 


E{P'} 


1 

2 

E-  , .  2n.‘ 

nr  (1-slno  sr> 

2 

2 

(81) 


E/No  is  the  expected  value  of  signal  energy  to  noise  power 
density  ratio.  The  predicted  normalized  mean-squared  error 
for  linear  phase  transition  and  Rayleigh  amplitude  is 


MSEN  =  §2  (8  - 


§> 


2 


1 

(1-sinc 

2 


2ir  '  I  2 


(82) 


The  results  of  equation  (82)  are  shown  for  various  values 
of  M  and  E/N0  in  Figures  (16)  and  (17).  In  general,  the 
predicted  MSE  for  the  Rayleigh  amplitude  case  is  higher  than 
the  MSE  for  the  corresponding  constant  amplitude  case.  As  in 
the  constant  amplitude  cases,  the  MSEN  becomes  smaller  as  the 
signal  to  noise  (E/N0)  is  increased  and  as  the  number  of  quan¬ 
tization  levels  is  increased. 

Step  Phase  Transition.  For  the  step  transition  case 
g(t)=6(t)  .  When  the  error  event  is  minimum  length  (N=2) 

and  ^ »  the  euclidean  distance  squared  is 


2n . 


D2  =  2E  i  -2j|iCOS(^-(a  i-8 1 )  ) 


(83) 


The  minimum  Dz  occurs  when  |cti-6i|  =  1 
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LINEAR  PHASE 


Figure  16  Predicted  MSEN  versus  M  ,  Linear  Transition,  Random  Variable 
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Figure  17  Predicted  MSEN  versus  E/No  ,  Linear  Transition,  Random  Variable 


The  expected  value  of  probability  of  any  error,  E^.,}  , 

£ 


is  given  by 


± 

E{p£}  =  E{P'}  =  E  )  fte  4N0 


m 


(84) 


e{pe}  = 


2+2  ^  d-cos  jf) 


(85) 


The  predicted  normalized  MSE  for  the  step  transition  case 
is  given  by 


MSEN 


1 

M2 


It  <8" 


2+2 


E 

No 


2tt  n 

(1-COS  jp) 


(86) 


The  results  for  this  case  are  shown  for  various  values  of  M 
and  E/N0  in  Figures  (18)  and  (19).  The  familiar  trend  to¬ 
wards  smaller  MSEN  with  increasing  E/N0  and  M  is  evident 
in  these  figures.  Overall,  the  predicted  MSEN  is  better  for 
this  case  than  for  the  linear  phase  transition  case.  The 
reason  for  this  is  that  there  is  a  larger  euclidean  distance 
between  two  signals  which  make  step  transitions  than  two  sig¬ 
nals  which  make  linear  transitions  (over  the  same  sequence  of 
phase  values). 

Bias 


It  has  been  shown  by  Scharf  (Ref  9:21-23)  that  the  max¬ 
imum  likelihood  (ML)  estimator  is  modulo-2iT  unbiased  for  the 
complex  normal  model  when  the  phase  is  constant.  As  pointed 


Figure  18  Predicted  MSEN  versus  M,  Step  Transition,  Random  Variable 
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Predicted  MSEN  versus  E/No,  Step  Transition,  Random  Variable 
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out  in  the  reference,  this  is  identical  to  the  MAP  case  with 
uniform  prior  statistics  on  the  phase  process. 

Performance  Analysis  Comments 

The  lower  bounds  in  figures  (9),  (11),  (13)  and  (15)  are 
strict;  however,  the  tightness  of  the  bounds  and  the  accuracy 
of  the  predicted  MSEN  depends  on  the  argument  of  the  Q  func¬ 
tion  in  the  prediction  equation.  Constant  argument  lines  are 
plotted  on  these  graphs.  There  is  a  high  degree  of  confidence 
in  the  predicted  values  whenever  the  argument  is  larger  than 
about  three  because  equation  (48)  is  a  good  approximation  to 
P£  in  this  region. 

As  can  be  seen  from  the  figures,  many  of  the  predicted 
points  are  outside  of  the  high  confidence  region.  Although 
the  predictions  may  be  good,  there  is  no  way  to  estimate  how 
good  they  are  with  these  methods.  In  addition,  these  perform¬ 
ance  predictions  give  no  indication  of  the  sensitivity  of  the 
estimator  to  decision  depth,  phase  variance,  or  various  corre¬ 
lation  times. 

Because  of  these  limitations  in  the  performance  analysis, 
a  simulation  of  the  estimator  is  useful. 


< 
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IV.  Estimator  Simulation 

A  Monte  Carlo  simulation  of  the  Viterbi  algorithm  modulo- 
2ir  phase  sequence  estimator  is  useful  to  obtain  performance 
data  for  an  entire  range  of  input  data  and  estimator  parameters. 

Simulation  Setup 

The  simulation  program  was  written  in  FORTRAN  and  run  on 
the  Cyber  175/Cyber  74  systems  available  for  use  at  The  Air 
Force  Institute  of  Technology.  After  the  parameters  for  the 
simulation  have  been  input,  the  program  generates  a  true  phase 
sequence,  constructs  the  signal  model  in  quadrature  form,  adds 
white  gaussian  noise  to  the  signal,  and  then  estimates  the 
phase  based  on  the  noise  corrupted  measurement.  The  phase  se¬ 
quence  is  estimated  over  one  scan  (100  measurements)  and  the 
error  and  mean-square  error  are  calculated  for  each  time  inter¬ 
val.  The  value  of  100  measurements  per  scan  was  selected  to 
be  much  larger  than  the  largest  decision  depth  (10)  so  that 
any  transient  effects  due  to  startup  could  be  observed.  Es¬ 
timations  can  be  repeated  for  many  runs  (scans)  in  order  to 
calculate  ensemble  averages  of  error  and  mean- square  error  at 
each  time  interval. 

Two  cases  of  true  phase  of  the  signal  are  used  in  the  sim¬ 
ulation.  In  the  first,  the  phase  is  constant  over  the  length 
of  the  scan  and  the  value  for  each  run  is  selected  randomly 
from  a  uniform  distribution  on  ( — tt  , tv )  .  This  allows 
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exercising  of  the  estimator's  parameters  without  any  masking 
effects  due  to  changing  phase.  The  random  choice  of  phase  is 
made  with  the  use  of  the  intrinsic  and  subroutine  functions 
for  selecting  uniform  random  deviates. 

The  second  sample  phase  function,  a  tactical  phase,  is 
used  to  approximate  what  could  occur  in  actual  use.  It  is 
assumed  that  the  dwell  time  T  is  0.5  ps  ,  the  illuminated 
area  is  a  circle  one  foot  in  diameter,  and  the  range  ambiguity 
distance  is  20  feet.  The  phase  of  a  signal  received  from  scan¬ 
ning  over  the  length  of  a  vehicle  (an  interesting  tactical  tar¬ 
get)  on  flat  ground  can  be  approximated  by  the  phase  function 
shown  in  Figure  20. 

The  phase  in  Figure  20  is  shown  beginning  at  -tt/2  radians 
for  clarity,  but  the  phase  is  allowed  to  start  at  any  value  sel¬ 
ected  from  a  uniform  distribution  on  (-1T,1T)  .  The  sample 

phase  in  Figure  20  is  modified  by  adding  zero  mean  white  gaus- 
sian  noise  with  a  variance  that  corresponds  to  one  foot  squared. 
This  is  done  to  depict  the  roughness  of  a  typical  surface  which 
would  be  found  under  real  conditions.  Instead  of  white  noise, 
the  phase  can  be  varied  by  adding  correlated  noise  in  order  to 
match  the  phase  statistics  of  a  random  walk  on  a  unit  circle. 

A  simple  way  to  do  this  is  by  filtering  white  gaussian  samples 
with  a  one  stage  infinite  impulse  response  (HR)  digital  fil¬ 
ter  and  then  adding  the  correlated  phase  noise  to  the  phase 
function  of  Figure  20. 

After  the  true  phase  sequence  has  been  constructed  for 
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the  whole  scan,  individual  measurements  are  calculated 

using  equation  (10).  Two  cases  of  amplitude  values  are  used 
in  the  results  that  follow.  In  the  first  case,  the  amplitude 
is  constant  over  the  entire  scan  and  the  value  of  A,  is  cal 

K 

culated  by 


.yfi 


£  2N0 
0 _ 


(87) 


In  the  second  case,  the  discrete  time  amplitude  values, 

,  are  selected  from  a  Rayleigh  distribution  to  simulate  a 
fading  signal.  The  second  moment  of  £  is  determined  by 
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where  E/N0  is  the  expected  value  of  the  signal  to  noise 
ratio.  The  second  moment  specifies  the  Rayleigh  distribution. 
Sample  values  are  calculated  by  generating  a  random  uniform 
sample  and  then  converting  it  to  a  Rayleigh  distribution  sam¬ 
ple.  The  coherence  time  of  the  fading  amplitude  can  be  varied 
by  letting  each  Rayleigh  sample  be  valid  for  n  measurement 
samples;  the  "coherence  time"  for  amplitude  is  then  approxi¬ 
mately  tc=nT 

In  the  following  results,  the  word  depth  refers  to  the 
decision  depth  of  the  algorithm.  This  is  identical  to  K0  , 
the  fixed  lag  described  in  Chapter  II.  At  time  k  ,  an  estim¬ 
ate  of  the  phase  at  time  k-K0  is  produced  by  the  algorithm. 

The  number  of  discrete  phase  values  the  estimate  can  attain  is 
labelled  LEVELS=M  in  the  figures  that  follow,  where  M  is  odd. 

The  term  a02T  from  equation(14)  is  labelled  VARPH  in  the 
graphical  results.  This  is  a  measure  of  the  variance  of  the 
phase  process  which  the  algorithm  "assumes"  it  is  tracking. 

As  o02T  becomes  small,  the  transition  density  of  equation  (14) 
becomes  gaussian  and  the  distribution  on  becomes 

N((J>k_i, o02T)  .  If  o02T  is  too  small,  the  estimator  will  not 

follow  rapid  fluctuations  in  the  true  phase.  Conversely,  if 
o02T  is  too  large,  the  estimator  will  be  susceptible  to  track¬ 
ing  isolated  noise  bursts  even  when  the  phase  is  constant. 
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The  number  of  ambiguous  decisions  for  a  scan  is  labelled 
AMBIG  in  the  figures.  This  is  a  figure  of  merit  which  gives 

A 

the  number  of  possible  values  of  <j>.  „  when  the  value  of 

K— XV  q 

A 

<J>k_K  for  each  of  the  M  sequences  is  not  identical.  In  gen¬ 
eral,  the  values  stored  in  the  M  sequences  will  converge  to 
a  single  sequence  as  time  progresses.  The  decision  rule  im¬ 
plemented  in  this  simulation  selects  the  sequence  with  the 
largest  and  then  outputs  the  k-K0  value  of  $  from 

that  sequence.  The  ambiguity  counter  is  incremented  once  for 

-  .  A 

each  sequence  which  does  not  have  the  same  value  for  <(>.  v 

K— K  q 

as  the  MAP  sequence.  An  ambiguous  decision  count  of  zero 
does  not  guarantee  that  the  correct  sequence  was  es¬ 

timated,  only  that  the  MAP  sequence  given  the  measurement 
sequence  was  estimated. 

Simulation  Results 

The  results  of  the  estimator  simulations  are  discussed  in 
this  section  and  shown  in  the  figures  that  follow.  The  parameter 
values  and  conditions  are  shown  on  the  figures. 

Confidence  of  Sample  Statistics.  The  data  presented  in  the 
following  figures  is  a  result  of  taking  ensemble  averages  of  the 
mean-squared  error  over  50  Monte  Carlo  runs.  This  value  was 
selected  by  observing  how  the  ensemble  mean  square  error  changes 
with  increased  number  of  runs,  shown  in  Figure(21).  As  the  num¬ 
ber  of  runs  increases,  the  ensemble  average  mean-square  error 
tends  to  settle  down  to  a  specific  value.  The  number  50  was 
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Figure  21  Number  of  Runs 


selected  to  provide  good  results  at  a  reasonable  expenditure 
of  computer  time.  Similar  results  can  be  obtained  by  observ¬ 
ing  the  ensemble  average  of  the  phase  errors  as  the  number  of 
runs  is  increased. 

The  weak  law  of  large  numbers  can  be  used  to  calculate  the 
probability  that  the  ensemble  average  of  the  phase  estimates 
is  more  than  0.1  radians  away  from  the  ensemble  average  of 
the  true  phase.  For  a  probability  of  0.1  that  the  phase 
difference  is  larger  than  0.1  radians,  3290  runs  are  nec¬ 
essary.  A  similar  argument  using  the  central  limit  theorem 
gives  a  result  that  809  runs  are  required. 

The  data  calculated  from  these  two  limit  theorems  fit 
the  data  collected  from  simulations.  Since  MSE  is  the  prime 
measure  of  performance  and  it  settles  quickly  for  number  of 
runs  greater  than  50,  the  value  of  50  was  selected  for  use  in 
the  generation  of  the  following  data. 

Sample  Estimation.  A  sample  phase  estimation  for  one 
scan  is  shown  in  Figure  (22).  The  normalized  ensemble  mean- 
squared  error  for  50  simulation  runs  with  these  parameters 
is  shown  in  Figure  (23).  The  true  phase  plot  is  marked  with 
a  "((>"  every  ten  points.  As  can  be  seen,  the  estimator  tracks 
the  true  phase  very  well. 

Because  the  estimator  is  modulo-2Tr  ,  estimates  may  be  ob¬ 
tained  which  appear  poor,  but  which  are  actually  good.  An  ex¬ 
ample  of  this  is  shown  in  Figures  (24)  and  (25).  In  this  ex¬ 
ample,  the  true  phase  crosses  over  the  +tt  boundary  several 
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Mean-Squared  Error 


Figure  24  Sample  Estimation 


Figure  25  Ensemble  Mean-Squared  Error 


times  between  N=40  and  60  .  The  estimator  produces  a  value 

"close  modulo-2ir"  to  the  true  value,  but  because  the  true 
phase  is  out  of  the  estimator's  output  range,  the  phase  estim¬ 
ate  appears  to  vary  wildly.  If  the  estimated  values  near  -it 
on  Figure(24)  are  adjusted  up  by  2tt  ,  the  figure  appears  very 
much  like  Figure  (22).  This  is  also  shown  by  comparing  the  en¬ 
semble  mean-squared  error  plots  corresponding  to  the  two  exam¬ 
ples.  The  mean-squared  errors  for  both  are  nearly  identical 
despite  the  fluctuations  in  the  estimated  phase.  This  is  anal¬ 
ogous  to  cycle  slipping  in  a  phase  locked  loop  (PLL).  This 
characteristic  of  a  modulo-2iT  estimator  must  be  considered 
when  designing  the  equipment  which  will  use  the  phase  estimate. 

Ambiguous  Decisions.  The  number  of  ambiguous  decisions 
varied  from  1400  for  depth=0,  to  700-800  when  depth=5  ,  and 
zero  when  depth=10  .  The  number  of  ambiguities  drops  quickly 
to  zero  as  depth  is  increased.  Very  little  variation  in  the 
number  of  mbiguities  was  found  as  other  parameters  such  as 
M,E/N0 ,  or  a02T  were  changed.  With  a  depth  of  10,  E/No>2  > 
and  either  fixed  or  varying  amplitude;  there  were  no  ambiguous 
decisions  observed. 

Parameter  Variations.  The  performance  of  the  estimator 
as  the  key  parameters  are  varied  is  displayed  in  Figures  (26) 
through  (42).  Analysis  of  the  individual  figures  is  given  in 
the  following  paragraphs.  Performance  is  graphed  as  normal¬ 
ized  mean-squared  error  versus  F/N0  and  families  of  curves 
are  plotted  for  varying  levels,  a02T  ,  depth,  amplitude  co¬ 
herence  time,  and  errors  in  the  estimation  of  the  amplitude. 
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Curves  are  plotted  for  the  four  combinations  of  constant 
phase,  constant  amplitude;  constant  phase,  Rayleigh  ampli¬ 
tude;  tactical  phase,  constant  amplitude;  and  tactical  phase, 
Rayleigh  amplitude.  The  figures  are  labeled  with  the  parameters 
of  interest  and  the  performance  of  the  estimator  is  evident. 

Figures  (26)  through  (29)  are  of  particular  interest  be¬ 
cause  they  correspond  to  the  parameters  used  in  predicting 
the  performance  of  the  estimator.  In  these  figures,  change 
in  performance  is  displayed  as  the  number  of  quantization  le¬ 
vels  is  varied.  Values  of  predicted  MSEN  and  the  lower  bound 
on  MSEN  (from  Chapter  III)  are  plotted  on  Figures  (26)  through 
(29)  for  the  case  of  M=15  .  A  comparison  of  the  predicted 

values  with  the  measured  results  shows  that  the  prediction  is 
good  in  areas  of  the  curve  where  the  argument  of  the  Q( • ) 
function  is  larger  than  three.  Also,  the  lower  bound  is  fairly 
tight.  This  is  as  expected  for  the  step  phase  transition  case 
since  the  Viterbi  algorithm  estimator  implemented  in  this  sim¬ 
ulation  has  as  its  input  discrete  phase  values  which  are  equiv¬ 
alent  to  a  continuous  phase  input  with  step  changes  in  value. 

As  predicted,  the  mean-squared  error  decreases  as  the  number  of 
quantization  levels  is  increased.  The  match  between  predicted 
MSEN  and  the  simulation  results  is  least  good  for  high  values 
of  M  and  low  values  of  E/N0  .  This  is  not  unexpected,  since 
the  argument  of  the  Q  function  is  small  for  low  E/N0  and  large 
M  . 

As  predicted,  the  performance  of  the  estimator  degraded 
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VARYING  LEVELS 


Figure  26  MSEN,  Varying  Levels, 
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Figure  27  MSEN,  Varying  Levels 
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Figure  28  MSEN  versus  E/No,  Varying  Levels,  Tactical 


Figure  29  MSEN  versus  E/No,  Varying  levels,  tactical 


when  the  amplitude  was  allowed  to  vary.  Simulation  results 
for  the  Rayleigh  fading  amplitude  case  for  various  values  of 
M  are  shown  in  Figures  (27)  and  (29).  Results  of  several  sim¬ 
ulations  indicate  that  the  crossover  of  the  M=15  and  M=31 
lines  is  a  result  of  the  randomness  of  the  Monte  Carlo  sim¬ 
ulation  for  a  sampling  of  50  runs  and  not  an  inherent  trait 
of  the  estimator.  These  results  also  correspond  closely  with 
the  predicted  MSEN  shown  in  Figure  (19)  for  the  random  ampli¬ 
tude  case. 

Figures  (30)  through  (33)  reflect  the  variation  of  mean- 
squared  error  performance  for  several  values  of  a02T  with 
M=15  .  As  described  previously,  a02T  is  the  variance  the 
estimator  "assumes"  for  the  transition  probabilities  of  the 
phase  process.  Figures  (30)  and  (31)  reflect  the  performance 
when  the  phase  is  constant  but  unknown  over  the  whole  scan. 

The  best  performance  is  achieved  when  a02T  is  small 
(a02T=O.l  for  the  values  displayed).  This  is  an  intuitive¬ 
ly  pleasing  result  since  a  constant  phase  has  a  conditional 
transition  variance  of  zero  and  this  MAP  estimator  will  tend 
to  select  constant  values  of  $k  when  a02T  is  very  small. 

It  should  be  noticed  that  performance  varies  more  at  low  val¬ 
ues  of  E/N0  than  at  high  values.  This,  too,  is  as  expected 
because  when  the  signal  amplitude  is  very  small,  the  estimator 
"relies"  more  heavily  on  the  transition  probabilities  than  the 
signal  measurement  for  calculation  of  the  rk’s  . 

The  performance  is  reversed  for  cases  in  which  the  phase 


Figure  30  MSEN,  Varying  VARPH, 


Figure  31  MSEN, 


Figure  32  AISEN,  Varying  VARPH,  Tactical 
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Figure  33  MSEN,  Varying  VARPH 


is  allowed  to  vary  as  shown  in  Figures  (32)  and  (33).  In 
these  simulations,  white  noise  samples  with  a  variance  of 
0.31  radians  squared  were  added  to  the  underlying  tactical 
phase  shape.  As  can  be  seen  from  these  figures,  as  a02T 
was  increased  from  0.1  to  0.5  radians  squared,  the  MSEN  per¬ 
formance  improved.  Further  increase  in  a02T  produced  no 
significant  improvement  in  performance.  A  large  value  of 
o02T  may  allow  the  estimator  to  track  the  phase  of  a  noise 
burst  instead  of  the  true  phase  of  the  signal  for  a  specific 
realization  of  the  random  measurement  process;  however,  this 
degradation  in  performance  as  a02T  becomes  too  large  is  not 
evident  in  the  figures  because  they  show  results  of  50  scans 
and  the  data  are  averaged  over  the  50  scans. 

Figures  (34)  through  (37)  display  the  performance  of  the 
estimator  as  the  decision  depth  is  varied.  In  all  of  these 
simulations,  the  number  of  quantization  levels  is  15.  The 
results  closely  correspond  to  the  results  predicted  in  Chapter 
III.  There  is  little  variation  in  performance  due  to  a  change 
in  depth  over  the  domain  of  signal  to  noise  ratios  displayed. 
Figure  (35)  does  show  a  tendancy  of  lower  MSEN  for  larger  val¬ 
ues  of  depth. 

Figures  (38)  and  (39)  show  the  effect  of  varying  the  coher¬ 
ence  time  of  the  random  amplitudes.  As  defined  previously,  "co¬ 
herence  time"  is  t  =nT  ,  where  T  is  the  sampling  time  and 
n  is  the  number  of  time  intervals  for  which  a  specific  realiza¬ 
tion  of  the  random  amplitude  is  valid.  There  is  little  change 
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Figure  34  MSEN,  Varying  Depth, 
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Varying  Depth,  A,  Tactical 


Figure  37  MSEN,  Varying  Depth,  A,  Tactical 


Figure  38  MSEN,  Varying  Coherence  Time, 


VARYING  COH.  TINE 


The  per- 


ff 


in  performance  as  t  varies  from  T  to  10T . 

G 

for^kance  curves  match  closely  the  predicted  performance. 

Figures  (40)  and  (41)  display  the  effects  of  perturba¬ 
tions  in  the  estimation  of  the  discrete  time  signal  amplitudes, 
.  The  perturbed  A£  is  calculated  for  the  inphase  case  by 

+  <89> 

and  for  the  quadrature  case  by 


&  =  V(Ak+fci>2+* 


(90) 


where  Afc  is  the  Rayleigh  amplitude,  and  and  a^  are 

Gaussian  random  variables  with  zero  mean  and  variance  a  T2 

a  I 

and  ,  respectively.  Little  change  in  performance  is 

exhibited  by  the  estimator  over  the  range  of  disturbed  ampli¬ 
tude  samples.  The  estimator  MSEN  increases  by  a  factor  of 
two  as  the  variance  of  the  disturbance  is  increased  from  0.1 
to  2.0  (at  E/No=2),  The  change  in  MSEN  is  less  at  higher 
values  of  E/N0  .  In  general,  the  estimator  is  insensitive 
to  disturbances  in  the  amplitude  estimation.  This  is  an  ex¬ 
tremely  useful  result.  Not  only  is  the  estimator  insensitive 


to  a  fast  fluctuating  signal,  it  is  not  greatly  sensitive  to 
errors  in  the  amplitude  estimate. 

Figure  (42)  shows  the  results  of  the  estimator  performance 
when  the  input  is  as  shown  in  Figure  (20)  with  correlated  noise 
added.  The  "noise"  portion  of  the  input  phase  was  calculated 
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Figure  40  MSEN  versus  E/N0  ,  In  Phase  Amplitude  Estimation  Disturbance 
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Figure  41  MSEN  versus  E/N0  ,  Quadrature  Amplitude  Estimation  Disturbance 


ni-m , 


recursively  by 


*k  "  *k-i  +  4  <91> 

where  is  a  Gaussian  random  variable  with  zero  mean  and  var¬ 
iance  of  0.5 

The  results  shown  in  Figure  (42)  are  very  close  to  those 

of  Figure  (33)  (same  parameters  except  the  input  phase  is 

"roughened"  with  white  noise).  The  estimator  displayed  little 

2 

sensitivity  to  the  match  between  the  ao  T  of  the  transition 
density  and  the  variance  of  the  input  phase  process. 

The  simulation  results  compare  favorably  with  the  per¬ 
formance  predictions  developed  in  Chapter  III.  In  addition,  the 
simulation  results  provide  information  on  the  effect  of  the 
paramenuers  depth,  a*T  ,  coherence  time,  and  amplitude  dis¬ 
turbance  on  performance.  The  next  chapter  compares  these  re¬ 
sults  with  other  estimation  schemes. 


I 
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V.  Comparison  with  Other  Estimators 

A  comparison  of  the  performance  of  this  estimator  with 
bounds  on  estimator  performance  and  predicted  results 
of  other  methods  is  necessary  in  order  to  judge  which  scheme 
is  best. 

Cramer  Rao  Lower  Bound 

Rife  and  Boorstyn  (Ref  8)  have  shown  that  the  variance  of 
a  phase  estimate,  $  ,  can  be  lower  bounded  by  the  Cramer  Rao 

(CR)  bound 

2 

Var(4>}  >  -2 —  (92) 

bfl2N 

where  b0  is  the  amplitude  of  the  real  signal  b  ocos(  wt+  <J>)  , 

a2  is  the  variance  of  additive  zero  mean  Gaussian  noise,  and 
N  is  the  number  of  discrete  time  observations  of  the  same 
constant  phase  signal  made  by  the  estimator.  This  CR  bound 
is  for  a  maximum  likelihood  (ML)  estimator  and  corresponds  to 
the  Viterbi  algorithm  estimator  when  the  a  priori  statistics 
on  phase  are  uniform  over  the  range  -it  to  it  .  The  amplitude 
(b0)  and  radian  frequency  (w)  are  assumed  constant  and  known. 

A 

When  N=1  ,  one  observation  is  made  for  each  value  of  <J> 

estimated.  This  corresponds  to  a  decision  depth  of  zero  in 
the  Viterbi  algorithm  estimator. 

For  the  signal  model  posed  in  Chapter  II,  the  signal  to 
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noise  ratio  of  the  input  signal  is 


E  _  AK2T 
No  2N  o 

and  the  variance  of  the  additive  noise  is 


(93) 


The  CR  bound  (normalized  by  tt2/3)  is 


(94) 


Var(4>} 


(95) 


This  CR  lower  bound  is  shown  in  Figure  (43)  along  with 
simulation  results  of  the  ML  Viterbi  algorithm  estimator  for 
constant  amplitude  and  phase  with  depth  equal  to  zero.  This 

A 

assumes  that  the  estimator  is  unbiased  so  that  Var{<}>}  =  MSE. 

The  CR  lower  bound  for  a  larger  N  can  be  easily  calcul¬ 
ated  from  equation  (87);  however,  the  problem  posed  by  Rife 
and  Boorstyn  is  that  of  estimating  constant  but  unknown  par¬ 
ameters.  The  estimation  problem  investigated  in  this  thesis 
is  that  of  phase  sequence  estimation  in  which  the  phase  can 
assume  different  values  at  each  discrete  time.  Thus,  the  low¬ 
er  bound  given  by  equation  (92)  for  N>1  is  not  relatable  to 
results  obtained  from  a  simulation  of  the  estimator  presented 
in  this  report. 
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Figure  43  CR  Lower  Bound  and  ML  Estimator  Performance 


Phase -Locked-Loop,  Constant  Amplitude 

The  performance  of  a  phase -locked- loop  for  phase  demodu¬ 
lation  is  well  modeled  and  useful  for  comparison. 

The  phase -locked-loop  (PLL)  is  assumed  to  be  first  order 
and  linearized.  The  variance  of  the  phase  estimate  is  given 
by 


where  E/N0;L  is  the  input  signal  to  noise  ratio  above  threshold. 

Values  of  o2^  ,  normalized  by  tt 2 / 3  are  shown  in 

Figure  (44)  plotted  with  simulation  results  for  various  numbers 
of  quantization  levels.  The  PLL  results  have  been  taken  from 
Viterbi  (Refl2:94). 

When  compared  with  the  simulation  case  of  15  quantization 
levels,  the  following  results  can  be  observed.  The  Viterbi  es¬ 
timator  outperforms  the  PLL  by  almost  six  db  at  low  signal  to 
noise  ratios  and  by  about  two  db  at  the  higher  signal  to  noise 
ratios  displayed  in  this  :igure. 

It  is  also  evident  from  Figure  (44)  that  the  PLL  outper¬ 
forms  the  M=3  and  M=7  Viterbi  estimator  above  a  specific 
"crossover"  signal  to  noise  ratio.  This  is  due  to  the  quanti¬ 
zation  MSE  lower  bound  inherent  in  the  discrete  value  represent¬ 
ation  of  the  phase  by  the  Viterbi  estimator.  The  other  two 
cases  (M=15  and  M=31)  also  have  a  lower  bound  on  MSE;  however, 
the  bound  is  reached  at  a  higher  signal  to  noise  ratio  than 
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Figure  44  Simulated  Estimator  Performance  versus  PLL,  Constant 


displayed  in  this  figure. 


Phase-Locked-Loop ,  Rayleigh  Amplitude 

The  predicted  results  for  PLL  performance  have  been  ex¬ 
tended  by  Weber  (Ref: 13)  to  include  the  case  of  phase  extima- 
tion  for  a  signal  with  Rayleigh  slowly  fading  amplitude. 

Results  for  the  slowly  fading  amplitude  case  are  shown  on 
Figure  (45).  All  conditions  are  the  same  as  for  the  PLL  con¬ 
stant  amplitude  case  with  the  exception  that  the  amplitude  is 
now  allowed  to  fade. 

A  comparison  with  the  M=15  case  can  again  be  made  show¬ 
ing  that  the  PLL  requires  an  input  signal  t’o  noise  ratio  of 
approximately  20  (13db)  to  achieve  the  same  mean-squared  error 
the  Viterbi  estimator  attains  at  a  signal  to  noise  ratio  of  two 
(three  db).  This  10  db  margin  is  nearly  constant  over  the  do¬ 
main  of  signal  to  noise  ratios  displayed  in  Figure  (45).  Weber 
also  notes  that  as  the  fading  bandwidth  becomes  large  (or  the 
loop  bandwidth  bandwidth  becomes  small)  the  phase  tracking  be¬ 
comes  worse  until  the  PLL  loses  lock  (Refl3:497). 

Comments 

From  the  foregoing  comparisons,  it  can  be  seen  that  in  the 
constant  amplitude  case  with  E/N0>lO  ,  the  Viterbi  estimator 
(M=15)  outperforms  the  PLL  by  up  to  5  db.  This  margin  decreases 
with  increasing  E/N  until  the  Viterbi  estimator's  performance  is 
limited  by  quantization  noise.  Below  a  signal  to  noise  ratio  of 
10,  the  smoothing  effect  and  use  of  prior  in  the  MAP  Viterbi 
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Figure  45  Simulated  Estimator  Performance  versus  PLL,  Slowly  Fading  A 


estimator  provide  significant  improvements  over  the  PLL. 

Similar  results  are  evident  for  the  fading  amplitude 
case  when  the  quantization  is  over  15  levels.  The  performance 
margin  of  the  Viterbi  estimator  is  approximately  10  db  better 
than  the  PLL.  It  is  important  to  note  that  the  Viterbi  algo¬ 
rithm  MAP  estimator’s  performance  is  achieved  under  conditions 
in  which  the  amplitude  is  allowed  to  vary  from  sample  to  sam¬ 
ple.  The  Rayleigh  fading  PLL  results  were  limited  to  a  slow¬ 
ly  fading  signal.  It  can  be  expected  that  the  PLL  would  per¬ 
form  much  worse  under  the  rapid  fading  conditions  and  eventually 
lose  lock  as  the  fading  bandwidth  is  increased. 


VI .  Implementation  Considerations 

This  chapter  will  consider  the  storage  requirements,  com¬ 
plexity,  and  speed  requirements  for  implementation  of  the  Vit- 
erbi  algorithm  for  modulo-2  ir  phase  sequence  estimation.  The 
problem  of  correct  interpretation  of  modulo-2ir  data  will  also 
be  discussed. 

Storage  Requirements 

K 

The  estimator  must  store  M  sequences,  {£k}i  ,  one  for 

each  of  the  M  phase  states.  Each  sequence  contains  Ko  +1 
floating  point  numbers  representing  the  estimates  of  phase  over 
the  time  interval  (k-Ko  )T  to  kT  .  This  is  a  total  of 
(M) • (Kq  +1)  locations.  For  each  sequence  a  value  of  r v  must 
be  stored.  The  current  time  counter,  k  ,  must  be  stored  as 
well  as  the  values  of  conditional  transition  probabili¬ 
ties.  This  gives  a  total  of  Mko  permanent  storage 

locations.  Temporary  storage  is  necessary  for  the  estimated 
signal  amplitude  ,  the  measurement  amplitude,  ,  and 

the  measurement  argument ,  <j>  k  ,  Temporary  storage  will  also 

be  required  for  holding  the  values  of  T  during  calculation 

Km 

of  the  algorithm.  This  requirement  can  vary  greatly.  For  ex¬ 
ample,  if  all  the  M  values  of  are  computed  in  parallel, 

m 

then  all  M  values  must  be  stored  until  the  biggest  is  picked. 

On  the  other  hand,  if  the  values  of  are  computed  serial- 

m 

ly,  then  only  the  current  and  "biggest-to-date”  values  need  be 
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stored.  Because  of  the  large  time  versus  storage  trade  off 
possible  with  temporary  storage,  more  definite  statements  on 
the  requirements  are  not  possible  until  the  specific  imple¬ 
mentation  of  the  algorithm  is  selected. 

Complexity 

This  complexity  and  speed  analysis  will  concern  only  the 

serial  algorithm  stated  in  Chapter  II.  It  is  assumed  that  a 

separate  processing  circuit  can  select  the  output  value  <j>.  v 

k-K 

independent  of  the  algorithm  computation.  It  is  also  assumed 
that  all  possible  calculations  are  made  outside  of  the  repet¬ 
itive  algorithm  to  conserve  time. 

In  order  to  advance  all  M  sequences  one  step  in  time,  M 
values  of  must  be  computed.  For  each  one,  M  values  of 

rk  must  be  computed  and  M-l  magnitude  comparisons  must  be 
m 

made.  The  calculation  of  one  value  of  rk  requires  two 

m 

additions,  one  subtraction,  one  cosine  calculation,  and  three 

multiplies.  Memory  reference  times  for  fetching  values  are 

ignored  since  dedicated  hardware  designed  for  the  Viterbi 

algorithm  can  greatly  optimize  these  times. 

If  the  times  required  for  the  above  calculations  are 

defined  as  follows: 

t„  *=  addition  or  subtraction  time 

a 

tm  *  multiply  time 

t _  «  comparison  time  (97) 

comp 

t„ _ ■  cos  (•)  calculation  time 


t 


then  the  total  time  to  calculate  one  value  of  is 

m 


t„  =  3t  +t  +3t 
r  a  cos  m 


(98) 


and  the  total  time  required  to  advance  all  M  sequences  by  one 
value  is 


l 


i 


< 


Li 


t„.  =  M2t_+M(M-l)t 
VA  r  comp 

(99) 

=  M2  (tr,  +  t  ) 

v  r  comp'' 

Therefore,  the  complexity  of  the  algorithm  is  KM2 

This  can  be  computationally  unfeasible  for  large  values  of  M, 

K 

but  is  certainly  better  than  the  M  complexity  previously 
described  for  the  global  search  version  of  MAP  sequence  estim¬ 
ation. 

Modulo-2iT  Data  Interpretation 

As  previously  described  in  the  results  section  of  Chapter 
IV,  a  phase  sequence  which  is  estimated  very  well  modulo-2ir 
may  appear  to  be  quite  bad. 

One  possible  solution  is  to  use  a  mapping  function  on  the 
estimated  phase,  $  ,  to  produce  the  display  phase  <J>D  .  Two 

example  functions  are  shown  in  Figure  (46). 

Both  mapping  functions  will  satisfy  the  requirement  that  a 
wraparound  of  <|>  from  it  to  -it  produces  very  little  change 
in  the  displayed  phase  sequence.  However,  both  methods  reduce 
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Fig.  46.  Display  Mapping  Functions 


the  range  ambiguity  by  a  factor  of  two  and  they  tend  to  accent 
uate  any  noise  in  the  phase  estimate.  In  addition,  the  cosine 
mapping  is  nonlinear. 

Another  possible  method  to  improve  the  display  character¬ 
istics  is  to  adjust  the  phase  estimate  up  or  down  by  2v  rad¬ 
ians  depending  on  certain  conditions.  One  method  for  doing 
this  can  be  described  by  the  algorithm: 


For  each  ^ ,  k=l,2,....,k 

Begin 

If  |||J  >  threshold  and  l^+^l  >  threshold 

Begin 

If  IV^k+i1  >  * 

Begi  n 

Adjust  $|,,  +  1  up  or  down  by  2ir  to  make  the 
the  sign  of  $^+1  equal  the  sign  of 

End 

End 

End 

This  algorithm  simply  adjusts  the  next  value  of  $  up 
or  down  by  2n  when  both  the  current  value  and  the  next  value 
are  close  to  the  tt  or  -tt  boundary.  As  the  threshold  ap¬ 
proaches  t r  ,  fewer  adjustments  are  made.  There  will  be 
cases  for  which  the  algorithm  is  "fooled"  into  making  an  un¬ 
warranted  adjustment,  but  these  occurrences  can  be  minimized 
by  proper  selection  of  the  threshold. 
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VII .  Conclusions  and  Recommendations 

This  report  has  dealt  with  the  problem  of  nonlinear  phase 
sequence  estimation  in  the  presence  of  additive  white  Gaussian 
noise  and  Rayleigh  fading.  Perfect  knowledge  of  each  realiza¬ 
tion  of  the  signal  amplitude  has  been  assumed.  The  problem 
was  motivated  by  the  laser  line-scan  imagery  system  develop¬ 
ed  by  the  Air  Force  Avionics  Laboratory;  however,  the  estima¬ 
tor  and  performance  indicators  derived  have  a  much  broader  ap¬ 
plication.  This  estimator  can  provide  excellent  phase  se¬ 
quence  estimates  for  other  communications  and  optics  problems 
in  which  the  phase  is  dynamically  changing  and  the  signal  am¬ 
plitude  is  rapidly  varying. 

Conclusions 

The  modulo-2TT  phase  sequence  Viterbi  algorithm  estimator 
with  15  phase  quantization  levels  outperforms  a  first  order  PLL 
operating  on  a  constant  amplitude  signal  by  2-6  dB  over  the  in¬ 
put  signal  to  noise  range  of  3-17  dB.  For  each  value  of  the 
number  of  quantization  levels,  there  is  a  value  of  signal  to 
noise  ratio  above  which  the  PLL  performance  is  better.  This 
is  due  to  the  inherent  quantization  noise  of  the  Viterbi  MAP 
estimator. 

For  rapidly  fading  signals,  the  Viterbi  MAP  estimator 
significantly  outperforms  a  first  order  PLL.  The  estimator 
requires  approximately  10  dB  less  input  signal  to  noise 


r 


ratio  to  equal  the  performance  of  a  PLL  estimating  the  phase 
of  a  signal  with  slowly  fading  amplitude.  As  the  fading  band¬ 
width  is  increased  for  the  PLL  estimator,  performance  rapidly 
degrades  until  the  PLL  loses  lock. 

Imperfect  knowledge  of  the  signal  amplitude  does  degrade 
the  phase  estimate  by  a  factor  of  about  two  when  E/N0=2  ; 
however,  there  is  little  degradation  at  higher  signal  to  noise 
ratios.  The  fact  that  the  estimator  is  insensitive  to  errors 
in  the  amplitude  estimate  is  significant  because  this  implies 
that  the  assumption  of  perfect  knowledge  of  the  signal  ampli¬ 
tude  is  not  critical  to  the  performance  predicted  in  this  re¬ 
port. 

The  performance  indicators  developed  in  Chapter  IV  and 
extended  to  the  Rayleigh  fading  case  have  been  validated  by 
the  Monte  Carlo  simulation  of  the  Viterbi  estimator.  The  pre¬ 
dicted  MSE  is  accurate  when  the  signal  to  noise  ratio  is  large 
or  when  the  number  of  phase  quantization  levels  is  small.  As 
shown  by  the  simulation  results,  the  predicted  results  are  al¬ 
so  useful,  although  not  as  accurate,  for  other  regions  of  the 
data. 

The  computational  complexity  of  the  algorithm  is  propor¬ 
tional  to  KM2  ,  where  M  is  the  number  of  discrete  phase 
levels  and  K  is  the  length  of  the  sequence.  This  compares 

ir 

favorably  with  the  M  complexity  of  a  global  search  for  the 
MAP  phase  sequence. 
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i  .  Recommendations  for  Further  Study 

The  following  areas  are  recommended  for  further  study. 

The  statistics  on  the  phase  process  only  assumed  a  corre¬ 
lation  in  the  direction  of  the  scan.  It  may  be  possible  to 
exploit  the  correlation  in  two  dimensions  to  improve  the  phase 
estimate. 

It  is  recommended  that  the  joint  estimation  of  ampli¬ 
tude  and  phase  be  investigated.  Solution  of  the  joint  prob¬ 
lem  will  eliminate  the  need  for  a  separate  amplitude  required 
by  the  Viterbi  estimator  described  in  this  report.  In  addi¬ 
tion,  it  may  be  possible  to  exploit  any  correlation  between 
amplitude  and  phase  in  order  to  improve  the  phase  estimate. 

» 

'  "  Filters  with  sharp  cutoffs  are  often  required  to  reduce 

the  aliasing  effects  described  in  Chapter  II.  These  filters 
often  have  nonlinear  phase  responses  as  a  result  of  the  re¬ 
quirement  for  a  sharp  cutoff.  It  is  recommended  that  a  model 
be  developed  which  includes  the  effects  of  aliasing.  With 
this  model,  it  can  be  determined  if  the  filtering  require¬ 
ments  can  be  made  less  stringent  while  still  maintaining  an 
acceptable  level  of  performance  for  the  phase  sequence  estim¬ 
ator. 

This  report  considered  phase  estimation  of  a  signal  in 
white  gaussian  noise.  In  optical  systems  such  as  this  lin- 
scan  imager,  quantum  noise  with  point  process  statistics  is 
a  significant  part  of  the  overall  system  noise.  It  is 
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recommended  that  a  similar  estimator  be  developed  for  a  signal 
which  includes  quantum  noise. 

The  simulation  of  this  estimator  has  assumed  either  per¬ 
fect  knowledge  of  the  signal  amplitude  or  perfect  knowledge 
corrupted  by  gaussian  quadrature  noise.  It  is  recommended 
that  the  actual  performance  be  determined  given  actual  ampli¬ 
tude  estimation  schemes. 
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In  this  report,  a  maximum-a  posteriori  (MAP)  phase  sequence 
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signal  is  observed  in  additive,  zero  mean,  white  Gaussian  noise 
and  the  statistics  of  the  amplitude  are  assumed  to  be  Rayleigh. 
Perfect  knowledge  of  each  realization  of  the  signal  amplitude 
is  assumed  and  the  sensitivity  to  errors  in  the  amplitude  estimate 
is  investigated.  Expressions  for  approximate  mean-squared  error 
and  a  lower  bound  on  mean-squared  error  are  developed  and 
verified  by  a  Monte  Carlo  simulation  of  the  estimator. 

For  quantization  of  the  phase  estimates  to  15  leveia  and 
for  a  signal  with  constant  amplitude,  the  estimator  outperforms 
a  first  order  phase-locked-loop  (PLL)  by  2-6  dB  over  the  range 
of  input  signal  to  noise  ratios  of  3-17  dB. 

For  the  case  when  the  amplitude  is  rapidly  fading,  the 
estimator  outperforms  a  first  order  PLL  by  greater  than  10  dB. 

This  performance  margin  is  based  on  the  Viterbi  estimator's 
performance  for  a  rapidly  fading  signal  compared  to  PLL 
performance  for  a  slowly  fading  signal.  Note  that  as  the  ampli¬ 
tude  fading  bandwidth  of  the  PLL  signal  is  increased,  the  PLL 
will  lose  lock  completely,  so  that  the  improvement  margin  present¬ 
ed  is  extremely  conservative. 

The  Viterbi  estimator  is  shown  to  be  insensitve  to  errors 
in  the  estimate  of  the  signal  amplitude  for  input  signal  to  noise 
ratios  larger  than  2  (  3  dB  ).  The  MAP  sequence  estimator 
is  also  shown  to  reduce  its  mean-squared  error  with  increasing 
decision  depth  (fixed  lag  of  the  Viterbi  algorithm),  and  increas¬ 
ing  number  of  discrete  phase  levels.  The  performance  depends 
weakly  on  the  match  between  the  variance  of  the  actual  phase 
and  the  variance  assumed  by  the  estimator. 
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